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ABSTRACT 


Unsteady  flow  past  three-dimensional  blunt  bodies  such  as  sphere,  rectangular  parallelepiped,  and 
circular  cylinder  has  been  numerically  analyzed  using  commercial  code,  Fluent,  to  offer  a  design 
guideline  for  Micro  Air  Vehicle.  Three-dimensional  Navier-Stokes  equations  are  solved  with  Large  Eddy 
Simulation  turbulence  model.  Three  different  aspect  ratios  of  1.0,  1.5,  and  2.0  are  selected  to  figure  out 
the  effect  of  geometric  shape  on  the  unsteady  flow  physics  in  the  range  of  Reynolds  numbers  from 
1.0  x  103  to  1.0  X  104.  The  mechanism  of  vortex  shedding  formation,  the  frequency  and  amplitude  of 
the  oscillating  unsteady  force  coefficients  including  drag,  lift  and  side  forces  are  investigated  through  this 
research.  There  are  some  differences  in  the  vortex  shedding  patterns  depending  on  the  blunt  body  shapes. 
The  frequency  and  the  amplitude  of  the  oscillating  unsteady  forces  strongly  depend  on  the  blunt  body 
shape,  aspect  ratio  of  the  blunt  body,  and  Reynolds  number. 


STATEMENT  OF  WORK 


The  objective  of  the  present  research  is  to  study  the  effects  of  three-dimensional  geometric  shapes 
and  Reynolds  numbers  on  the  unsteady  flow  physics  around  blunt  body,  which  is  flying  through  the  air, 
concentrating  on  the  vortex  shedding  mechanism,  oscillating  force  frequency  and  its  amplitude  exerted  on 
the  body  to  prepare  a  design  guideline  for  MAV.  This  kind  of  research  will  give  the  design  baseline  to 
cope  with  the  wind-induced  vibration  which  may  cause  out  of  control  or  destruction  of  MAV.  A  sphere,  a 
rectangular  parallelepiped,  and  a  circular  cylinder  are  selected  as  geometrically  simplified  MAVs  and 
numerical  analysis  of  the  unsteady  flowfield  around  these  blunt  bodies  have  been  accomplished  at  the 
Reynolds  numbers  of  1.0  X  103,  2.0  X  103,  4.0  X  103,  and  1.0  X  104.  Three  different  aspect  ratios  of 
1.0,  1.5,  and  2.0  are  considered  for  a  rectangular  parallelepiped  and  a  circular  cylinder. 
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I.  INTRODUCTION 

The  periodicity  of  the  wake  of  a  blunt  body  is  associated  with  the  formation  of  a  stable  street  of 
staggered  vortices.  Karman[l]  analyzed  the  stability  of  vortex  street  configurations  and  established  a 
theoretical  link  between  the  vortex  street  structure  and  the  drag  on  the  body.  Jordan  and  Fromm[2] 
investigated  oscillatory  drag,  lift,  and  torque  on  a  circular  cylinder  in  a  uniform  flow  at  Reynolds  numbers 
of  100,  400,  and  1,000  by  solving  vorticity-stream  function  formulation.  They  showed  the  dramatic  rise  of 
the  drag  force  coefficient  during  the  development  of  the  Karman  vortex  street.  A  detailed  study  of  the 
wake  structures  and  flow  dynamics  associated  with  two-dimensional  flows  past  a  circular  cylinder  is 
performed  by  Blackburn  and  Henderson[3]. 

It  would  be  very  valuable  attempts  to  study  the  flows  past  elliptic  cylinders  because  engineering 
applications  often  involve  flows  over  complex  bodies  like  wings,  submarines,  missiles,  and  rotor  blades, 
which  can  hardly  be  modeled  as  a  flow  over  a  circular  cylinder.  In  such  flows,  cylinder  thickness  and 
angle  of  attack  can  greatly  influence  the  nature  of  separation  and  the  wake  structure  [4]. 

In  1987,  Ota  et  al.[5]  investigated  a  flow  around  an  elliptic  cylinder  of  axis  ratio  1:3  in  the  critical 
Reynolds  number  regime,  which  extends  from  about  Re  =  8.5  X  104  to  3.12  X  105,  on  the  basis  of 
mean  static  pressure  measurements  along  the  cylinder  surface  and  of  hot-wire  velocity  measurements  in 
the  near  wake.  Nair  and  Sengupta[6]  solved  Navier-Stokes  equations  in  order  to  study  the  onset  of 
computed  asymmetry  around  elliptic  cylinders  at  a  Reynolds  number  of  1.0  X  104.  They  found  that  the 
ellipses  developed  asymmetry  much  earlier  than  the  circular  cylinder.  Kim  and  Sengupta[7]  studied 
unsteady  flow  past  an  elliptic  cylinder  whose  axis  ratios  are  0.6,  0.8,  1.0,  and  1.2  at  different  Reynolds 
numbers  of  200,  400,  and  1,000  to  investigate  the  unsteady  lift  and  drag  forces.  They  found  that  the 
elliptic  cylinder  thickness  and  Reynolds  number  could  affect  significantly  the  frequencies  of  the  force 
oscillations  as  well  as  the  mean  values  and  the  amplitudes  of  the  drag  and  lift  forces. 

Many  people  also  investigated  effect  of  incident  angles.  Patel[8]  investigated  the  incompressible 
viscous  flow  around  an  impulsively  started  elliptic  cylinder  at  0°,  30°,  45°  and  90°  incidences  and  at  the 
Reynolds  numbers  of  100  and  200.  Chou  and  Huang[9]  studied  the  unsteady  two-dimensional 
incompressible  flow  past  a  blunt  body  at  high  Reynolds  numbers  up  to  4.0  X  104.  They  considered  the 
aspect  ratio  and  angle  of  attack  as  controlled  parameters.  In  2001,  Badr  et  al.  [10]  used  a  series  truncation 
method  based  on  Fourier  series  to  reduce  the  Navier-Stokes  equations.  The  Reynolds  number  range  was 
up  to  5.0  x  103  and  axis  ratios  of  the  elliptic  cylinder  were  between  0.5  and  0.6,  and  angle  of  attack 
range  between  0°  and  90°.  They  showed  an  unusual  phenomenon  of  negative  lift  occurring  shortly  after 
the  start  of  motion.  Kim  and  Park[ll]  studied  numerically  to  figure  out  the  effects  of  elliptic  cylinder 
thickness  (thickness  to  chord  ratios  of  0.2,  0.4,  and  0.6),  angle  of  attack  (10°,  20°,  and  30°),  and  Reynolds 
number  (400  and  600)  on  the  unsteady  lift  and  drag  forces  exerted  on  the  elliptic  cylinder.  Through  this 
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study,  they  observed  that  the  elliptic  cylinder  thickness,  angle  of  attack,  and  Reynolds  number  are  very 
important  parameters  to  decide  the  unsteady  characteristics  of  the  lift  and  drag  forces. 

Since  many  engineering  structures  have  rectangular  shaped  cross  sections,  many  researchers  have 
carried  unsteady  flow  investigations  around  rectangular  cylinders  experimentally  or  numerically  in  view 
of  Reynolds  number,  angle  of  attack,  and  blockage  ratio  effect.  Davis  and  Moore[12]  solved  two- 
dimensional  incompressible  Navier-Stokes  equations  to  investigate  the  vortex  shedding  phenomenon 
around  a  rectangle  at  the  Reynolds  numbers  of  100  to  2,800  with  different  angle  of  attacks  and  rectangle 
dimensions.  They  found  that  the  properties  of  vortices,  lift,  drag,  and  Strouhal  number  are  strongly 
dependent  on  the  Reynolds  numbers.  Okajima[13]  investigated  experimentally  vortex  shedding 
frequencies  of  various  rectangular  cylinders  in  a  wind  tunnel  and  in  a  water  tank.  He  found  that  there  was 
a  certain  range  of  Reynolds  number  for  the  cylinders  with  the  width-to-height  ratios  of  two  and  three 
where  flow  pattern  abruptly  changed  with  a  sudden  discontinuity  in  Strouhal  number.  In  1993, 
Norberg[14]  measured  the  pressure  distributions  along  the  rectangular  cylinder  surface  at  angles  of  attack 
0°  to  90°  with  side  ratios  of  1.0,  1.62,  2.5,  and  3.0.  He  also  obtained  Strouhal  numbers  using  hot-wire  in 
the  near  wake  regions.  He  found  that  the  flow  showed  a  large  influence  of  both  angle  of  attack  and  side 
ratio  due  to  reattachment  and  shear  layer  /  edge  interactions.  Sohankar  et  al.[15]  calculated  unsteady  two- 
dimensional  flow  around  a  square  cylinder  at  incidences  between  0°  to  45°  and  Reynolds  numbers  of  45 
to  200.  They  used  an  SIMPLEC  algorithm  with  a  non-staggered  grid  arrangement.  They  found  that  the 
onset  of  vortex  shedding  occurred  within  the  interval  40  <  Re  <  55,  with  a  decrease  in  Reynolds  number 
with  increasing  angle  of  attack. 

Most  real  engineering  structures  have  three-dimensional  geometric  shapes  and  the  flow  around  a 
three-dimensional  bluff  body  is  of  great  interest  in  engineering  practice.  Therefore,  many  investigators 
have  also  studied  three-dimensional  flowfield  around  blunt  bodies  to  have  better  understanding  of  the 
flow  physics. 

Mittal [16]  simulated  unsteady  flow  around  a  sphere  to  observe  vortex  dynamics  in  the  wake  of  a 
sphere  in  the  Reynolds  number  range  from  350  to  650  using  Fourier-Chebyshev  spectral  collocation 
method.  He  found  that  the  sphere  wake  at  these  transitional  Reynolds  numbers  exhibits  multiple  dominant 
frequencies  and  the  non-linear  interaction  between  these  frequencies  leads  to  a  complex  evolution  of 
vortex  structures  in  the  near  wake.  Jones  and  Clarke[17]  simulated  flow  around  a  sphere  in  several 
different  flow  regimes  such  as  steady-state  laminar  flow  at  a  Reynolds  number  of  1.0  X  102,  unsteady 
laminar  flow  at  a  Reynolds  number  of  3.0  X  102,  turbulent  flow  with  laminar  boundary  layers  at  a 
Reynolds  number  of  1.0  X  104,  and  turbulent  flow  with  turbulent  boundary  layers  at  a  Reynolds  number 
of  1.0  X  106  by  using  Fluent  commercial  code.  Through  this  study,  it  is  found  that  Fluent  is  able  to 
accurately  simulate  the  flow  behavior  in  each  of  the  above  flow  regimes.  Johnson  and  Patel[18] 
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investigated  the  flow  past  a  sphere  numerically  and  experimentally  at  Reynolds  numbers  of  up  to  300. 
Steady  axisymmetric  flow  occurred  at  the  Reynolds  number  of  up  to  200.  For  Reynolds  numbers  of  210 
to  270,  a  steady  non-axisymmetric  flowfield  was  observed  whereas  unsteady  flowfields  were  observed  at 
the  Reynolds  numbers  greater  than  270.  At  the  Reynolds  number  of  300,  a  highly  organized  periodic  flow 
pattern  was  obtained  due  to  dominant  vortex  shedding  like  hairpin-shaped  vortices.  Tomboulides  et  al.[19] 
had  a  Direct  Numerical  Simulation  (DNS)  at  the  Reynolds  numbers  from  25  to  1,000  and  a  Large  Eddy 
Simulation  (LES)  at  20,000  for  flow  past  a  sphere.  They  found  two  early  bifurcations  of  the  sphere  wake, 
the  first  at  Reynolds  number  of  212  leading  to  a  three-dimensional  steady-state,  and  the  second  at 
Reynolds  numbers  between  250  and  285  resulting  in  a  unsteady  periodic  flow.  A  shear  layer  instability  is 
also  resolved  accurately  at  Reynolds  number  of  1,000.  Their  LES  showed  a  good  agreement  with 
experimental  results  both  in  terms  of  frequency  spectrum  and  drag  force  coefficient.  Tomboulides  and 
Orszag[20]  studied  transitions  that  occur  with  increasing  Reynolds  number  in  the  flow  past  a  sphere  by 
using  mixed  spectral  element  /  Fourier  spectral  method.  They  found  that  the  first  transition  of  the  flow 
past  a  sphere  is  a  linear  one  and  leads  to  a  three-dimensional  steady  flowfield  with  planar  symmetry.  The 
second  transition  leads  to  a  single  frequency  periodic  flow  with  vortex  shedding,  which  maintains  the 
planar  symmetry  observed  at  low  Reynolds  number.  As  the  Reynolds  number  increases  further,  the  planar 
symmetry  is  lost  and  flow  reaches  a  chaotic  state.  Constantinescu[21]  simulated  a  flowfield  around  a 
sphere  for  the  conditions  of  subcritical  and  supercritical  regimes  at  the  Reynolds  numbers  from 
1.0  X  104  to  1.0  X  106.  They  devoted  their  attention  to  assess  pressure  distribution,  skin  friction,  and 
drag  as  well  as  to  understand  the  vortex  dynamics  with  the  Reynolds  number  using  Detached  Eddy 
Simulation  (DES).  Taneda[22]  had  a  visual  observation  of  the  flow  past  a  sphere  at  Reynolds  numbers 
between  l.OxlO4  to  1.0X106  by  means  of  the  surface  oil  flow  method,  the  smoke  method  and  the 
tuft  grid  method  in  the  wind  tunnel.  They  found  that  the  wake  performs  a  progressive  wave  motion  at 
Reynolds  numbers  between  1.0  X  104  and  3.8  X  105  and  forms  a  pair  of  streamwise  line  vortices  at 
Reynolds  numbers  between  3.8  X  105  and  1.0  X  106 .  Sakamoto  and  Haniu[23]  experimentally 
investigated  vortex  shedding  from  spheres  at  Reynolds  numbers  from  3.0  X  102  to  4.0  X  104  in  a 
uniform  flow  using  hot-wire  technique  in  a  low  speed  wind  tunnel.  Also,  flow  visualization  was  carried 
out  in  a  water  channel.  They  classified  the  variation  of  the  Strouhal  number  with  the  Reynolds  number 
into  four  regions  and  found  that  the  higher  and  lower  frequency  modes  of  the  Strouhal  number  coexisted 
at  Reynolds  numbers  ranging  from  8.0  X  102  to  1.5  X  104.  Constantinescu  et  al. [24,25]  numerically 
simulated  the  subcritical  flow  at  a  Reynolds  number  of  1.0  X  104  over  a  sphere  to  compare  prediction  of 
some  of  the  main  physics  and  flow  parameters  from  solutions  of  the  unsteady  Reynolds-averaged  Navier- 
Stokes  equations  (URANS),  LES,  and  DES.  URANS  used  k  —  e,  k  —  a),  v2  —  f,  and  Spalart-Allmaras 
model  and  LES  used  dynamic  eddy  viscosity  model.  DES  is  a  hybrid  method  which  hires  Reynolds- 


averaged  Navier-Stokes  equations  (RANS)  near  the  wall  and  LES  in  the  wake.  DES  and  LES  showed 
better  agreement  with  measurements  than  URANS  although  all  of  the  techniques  showed  very  similar 
profiles  of  the  mean  velocity  and  turbulent  kinetic  energy  in  the  near  wake.  Sakamoto  and  Haniu[26] 
experimentally  investigated  the  formation  mechanism  and  frequency  of  vortex  shedding  from  a  sphere  in 
uniform  shear  flow  in  a  water  channel  using  flow  visualization  and  velocity  measurement  at  the  Reynolds 
numbers  of  2.0xl02  to  3.0xl03.  The  shear  parameter  defined  as  the  transverse  velocity  gradient  of 
the  shear  flow  was  varied  from  0  to  0.25.  They  found  that  the  critical  Reynolds  number  beyond  which 
vortex  shedding  occurred  was  lower  than  that  for  uniform  flow  and  decreased  approximately  linearly  with 
increasing  shear  parameter.  Also,  they  found  that  the  Strouhal  number  of  the  hairpin-shaped  vortex  loops 
became  larger  than  that  for  uniform  flow  and  increased  as  the  shear  parameter  increased.  Unlike  the 
detachment  point  of  vortex  loops  in  uniform  flow,  which  was  irregularly  located  along  the  circumference 
of  the  sphere,  the  detachment  point  in  shear  flow  was  always  on  the  high-velocity  side.  Achenbach[27] 
had  an  experimental  study  of  vortex  shedding  from  spheres  in  the  Reynolds  number  range  4.0  X  102  to 
5.0  X  106.  Strong  periodic  fluctuations  in  the  wake  were  observed  from  Reynolds  number  6.0  X  103  to 
3.0  X  105  whereas  periodic  vortex  shedding  could  not  be  detected  by  his  measurement  techniques 
beyond  the  Reynolds  number  of  3.7  X  105.  Taneda[28]  photographically  investigated  wakes  produced  by 
a  sphere  in  a  water  tank  at  Reynolds  numbers  from  5  to  300.  The  permanent  vortex -ring  began  to  form  in 
the  rear  of  a  sphere  at  Reynolds  number  of  24  and  to  oscillate  when  Reynolds  number  reached  about  130. 
Bakic  et  al.[29]  experimentally  investigated  turbulent  structures  of  flow  around  a  sphere.  The  mean 
velocity  field  and  turbulence  quantities  were  obtained  at  Reynolds  number  of  5.0  X  104  by  using  laser- 
Doppler  anemometry  in  a  small  low  speed  wind  tunnel.  Also,  flow  visualization  at  Reynolds  numbers 
between  2.2  X  104  and  4.0  X  105  were  performed  in  the  bigger  wind  tunnel  and  water  channel. 

Understanding  of  the  flow  physics  around  a  rectangular  parallelepiped  is  very  helpful  to  design 
buildings,  vehicles,  and  bridges,  etc.  Krajnovic  and  Davidson[30]  analyzed  the  flowfield  around  a  surface 
mounted  cube  at  the  Reynolds  number  of  4.0  X  104  using  LES  to  obtain  drag,  lift,  and  vortex  shedding 
frequency.  Also,  they  studied  coherent  structures  and  flow  features.  Krajnovic  and  Davidson[31]  had  a 
feasibility  study  of  use  of  LES  in  external  vehicle  aerodynamics.  LES  of  the  flow  around  simplified  car¬ 
like  shapes  at  the  Reynolds  number  of  4.0  X  104  could  give  the  knowledge  of  the  flow  around  a  car. 
They  simulated  the  flow  around  a  cube  and  the  other  of  the  flow  around  a  simplified  bus.  Iaccarino  and 
Durbin[32]  performed  unsteady  three-dimensional  RANS  simulations  with  v2  —  f  turbulence  model  to 
the  solution  of  the  flow  around  a  surface  mounted  cube  at  the  Reynolds  number  of  4.0  X  104.  The  flow 
around  a  cube  exhibited  a  strong  horseshoe  vortex  and  arch-shaped  vortex  in  the  near  wake.  Vengadesan 
and  Nakayama[33]  evaluated  three  different  sub-grid  scale  stresses  (SGS)  closure  LES  models  for 
turbulent  flow  over  a  square  cylinder  at  the  Reynolds  number  of  2.2  X  104.  They  were  conventional 
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Smagorinsky  model,  Dynamic  model,  and  one-equation  model.  It  was  concluded  that  a  one-equation 
model  for  subgrid  kinetic  energy  was  the  best  choice  from  the  viewpoint  of  affordable  computer  resources 
and  reasonable  turnaround  times.  Murakami  and  Mochida[34]  analyzed  unsteady  flowfields  past  a  two- 
dimensional  square  cylinder  at  the  Reynolds  number  of  2.2  X  104  and  compared  two-dimensional  and 
three-dimensional  LES  results.  Three-dimensional  LES  results  agreed  very  well  with  the  experimental 
results,  but  the  results  based  on  two-dimensional  computation  were  different  from  the  experimental 
results.  Also,  they  compared  various  turbulent  models  such  as  standard  k  —  £,  modified  k  —  £,  and 
Reynolds  stress  model.  Modified  k  —  £  turbulent  model  succeeded  in  reproducing  vortex  shedding  very 
well.  In  general,  LES  showed  the  best  agreement  with  the  experimental  results  although  it  took  a  great 
deal  of  CPU  time. 

Luo  et  al.[35]  simulated  particle-laden  wakes  of  a  circular  cylinder  with  Reynolds  number  ranging 
from  140  to  260  to  understand  three-dimensional  dispersion  of  particles  in  the  flow  around  a  bluff  body. 
They  developed  a  Lagrangian  tracking  solver  to  trace  the  trajectories  of  each  particle  in  the  non-uniform 
grid  system  and  observed  coherent  structure  and  vortex  dislocation  frequency.  Karlo  and  Tezduyar[36] 
performed  parallel  three-dimensional  computation  of  unsteady  flows  around  circular  cylinders  at 
Reynolds  numbers  of  300  and  800.  The  three-dimensional  features  were  weak  at  Reynolds  number  of  300 
whereas  three-dimensional  features  were  stronger  at  Reynolds  number  of  800.  Strong  three-dimensional 
features  arose  from  the  instability  of  the  columnar  vortices  forming  the  Karman  street.  They  also 
simulated  the  flowfield  at  the  Reynolds  number  of  1.0  X  104  with  LES  turbulence  model.  The  features 
were  very  similar  to  those  from  Reynolds  number  of  800,  however,  the  structures  were  more  diffuse  due 
to  the  increased  turbulence.  Zhao  et  al.[37]  studied  the  transition  of  the  flow  from  2-D  to  3-D  for  the  flow 
past  a  stationary  circular  cylinder  at  yaw  angles  in  the  range  of  0°  to  60°  at  the  Reynolds  number  of 
1.0  X  103  using  DNS  technique.  The  streamwise  vortices,  the  vortex  dislocation  and  the  instability  of  the 
shear  layer  were  observed  in  the  flow  visualization  as  well  as  by  the  numerical  analysis.  The  effects  of  the 
yaw  angle  on  wake  structures,  vortex  shedding  frequency  and  hydrodynamic  forces  were  investigated. 
Yeo  and  Jones [38]  investigated  the  three-dimensional  characteristics  of  the  fully  developed  flow  past  a 
yawed  and  inclined  circular  cylinder  at  the  Reynolds  number  of  1.4  X  105  using  DES.  Axial  lengths  of 
10,  20  and  30  times  of  circular  diameter  were  simulated.  They  found  that  the  swirling  flow  with  low 
pressure  on  the  cylinder  played  an  important  role  in  generating  multiple  moving  forces  and  the  flow  was 
highly  three-dimensional.  Persillon  and  Braza[39]  studied  the  transition  to  turbulence  of  the  flow  around  a 
circular  cylinder  by  a  three-dimensional  numerical  simulation  of  the  Navier-Stokes  equations  system  in 
the  Reynolds  number  range  100-300.  Wissink  and  Rodi[40]  performed  a  DNS  of  incompressible  flow 
around  a  circular  cylinder  at  Reynolds  number  of  3.3x10s.  A  significant  production  of  turbulence 
kinetic  energy  was  observed  inside  the  rolls  of  recirculating  flow  as  the  shear  layers  rolled-up. 
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Micro  Air  Vehicle  (MAV)  is  a  small  sized  light  autonomous  flying  machine  and  is  operating  at  low 
Reynolds  numbers  of  1.0  X  103  to  1.0  X  105.  At  such  Reynolds  number  ranges,  MAV’s  performance 
strongly  depends  on  the  laminar-to-turbulent  transition  and  unsteady  characteristics.  Many  unresolved 
research  areas  remain  in  the  low  Reynolds  number  aerodynamics  [4 1-42]. 

In  the  present  research,  three-dimensional  unsteady  flow  simulations  past  blunt  bodies  such  as  sphere, 
rectangular  parallelepiped,  and  circular  cylinder  have  been  performed  using  Fluent  which  is  very  well 
known  commercial  flow  solver  developed  by  Fluent,  Inc. 

Most  researchers  who  have  performed  three-dimensional  unsteady  flowfield  analysis  around  a 
rectangular  parallelepiped  or  a  circular  cylinder  have  blunt  bodied  which  are  fixed  on  the  plate  or  have 
long  span.  In  other  words,  they  have  three-dimensional  blunt  bodies  but  they  do  not  have  perfect  three- 
dimensional  flowfield.  The  present  research,  however,  analyzes  the  three-dimensional  flowfield  around  a 
blunt  body  which  is  flying  through  the  air  so  that  the  flowfield  is  thoroughly  three-dimensional  flowfield. 
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II.  FLOW  SOLVER 
II-l.  Solver  Overview 

In  this  research,  a  commercial  code  will  be  used  to  simulate  the  flowfield  around  a  blunt  body.  There 
are  many  commercial  CFD  codes  in  the  market  including  PHOENICS,  Fluent,  FLOW3D,  STAR-CD, 
ANSYS  CFX,  etc.  In  this  research,  Fluent  and  CFX  are  considered  as  candidate  CFD  tools  and  one  of 
them  will  be  selected  as  a  flow  solver  by  personal  decision  based  on  the  accuracy  and  computational 
efficiency. 

ANSYS  Fluent[43]  is  a  computational  fluid  dynamics  computer  code  developed  by  Fluent  Inc. 
Navier-Stokes  equations  are  solved  using  cell-centered  finite-volume  method.  Fluent  offers  several 
options  such  as  coupled  explicit,  coupled  implicit,  and  segregate  method  to  solve  the  governing  equations. 
The  grid  generation  module  called  Gambit  is  available  to  generate  two-dimensional  or  three-dimensional 
grid  structures  by  using  several  types  of  computational  cells  including  triangular,  quadrilateral, 
hexahedral,  tetrahedral,  pyramidal,  prismatic,  and  hybrid  meshes.  Although  Gambit  is  a  useful  tool  to 
generate  surface  meshes  as  well  as  volume  meshes,  it  is  preferred  to  use  Tgrid  to  make  a  three- 
dimensional  volume  meshes  because  it  can  handle  domains  with  a  large  number  of  computational  cells. 
Tgrid  is  also  developed  by  Fluent  Inc.  and  has  similar  interface  and  graphical  structure  to  that  of  Fluent. 
All  the  numerical  results  are  plotted  using  built-in  post  module  or  CFX-Post. 

ANSYS  CFX[44]  is  a  general  purpose  CFD  code,  which  combines  an  solver  with  pre  and  post 
processors.  CFX-Solver  solves  unsteady  Navier-Stokes  equations  in  conservation  form  using  coupled 
method.  It  uses  second  order  numerics  by  default,  ensuring  users  always  get  the  most  accurate  predictions 
possible.  CFX-Pre  can  import  mesh  files  produced  by  other  grid  generation  packages.  In  addition,  flow 
physics,  boundary  conditions,  initial  values  and  solver  parameters  are  specified  in  CFX-Pre.  Once  the 
solver  runs,  CFX-Post  generates  a  variety  of  graphical  objectives  through  interactive  post-processing  job. 

II-2.  Turbulence  Model  Overview 

Unsteady  flow  past  a  blunt  body  including  sphere,  rectangular  parallelepiped,  and  circular  cylinder  is 
turbulent  flow  in  the  Reynolds  number  range  of  1.0  x  103  to  1.0  x  104  [17,23,27,30,31,32,33,36,37,40]. 

Turbulence  is  time-dependent,  three-dimensional,  highly  non-linear  flow  phenomena.  Turbulence 
theory  states  that  the  eddies  vary  in  size;  larger  eddies  break  down  into  smaller  eddies.  This  process  of 
eddy  breakdown  transfers  kinetic  energy  from  the  mean  flow  to  progressively  smaller  scales  of  motion, 
which  is  known  as  the  energy  cascade.  At  the  smallest  scales  of  turbulent  motion,  the  kinetic  energy  is 
converted  to  the  heat  by  means  of  viscous  dissipation. 

The  simplest  and  most  straight  forward  way  to  describe  a  flow  is  to  solve  the  Navier-Stokes 
equations  directly  without  any  approximations  applied  in  the  calculation.  This  method  is  known  as  the 
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DNS.  DNS  resolves  all  the  significant  scales  of  turbulent  flow  down  to  the  Kolmogorov  scales,  which 
mean  scales  responsible  for  the  viscous  dissipation  of  energy  in  the  flow.  However,  DNS  requires  a  lot  of 
computation  time  as  well  as  huge  size  of  grid  especially  large  Reynolds  number  flows.  Thus  it  is  still  not 
practical  to  resolve  turbulent  flows  at  high  speed  flow  region  using  DNS  with  the  currently  available 
computer  capability. 

In  the  late  1980’s  and  up  to  the  beginning  of  1990’s,  RANS  type  modeling  is  mostly  used  to  simulate 
turbulent  flows  with  the  advent  of  increasing  power  of  digital  computers.  In  developing  the  governing 
equations  to  describe  turbulent  flow,  there  exist  fluctuations  in  the  flow.  In  the  Reynolds  Averaged 
approach  all  flow  variables  are  divided  into  mean  values  and  fluctuations  and  then  are  averaged  by  time 
so  that  fluctuating  components  are  removed.  In  the  Navier-Stokes  equations  the  time  averaging  yields 
new  terms  which  are  known  as  Reynolds  stresses  and  involve  mean  values  of  products  of  fluctuations. 
Unfortunately,  the  Reynolds  stresses  are  unknown  because  the  velocity  fluctuations  are  not  computed 
directly.  Therefore  turbulence  models  are  necessary  and  the  most  widely  used  turbulence  models  are 
Spalart-Allmaras  model[45],  k  —  a  model[46],  and  k  —  co  model.  During  the  last  few  years  URANS 
type  of  approach  has  been  also  attempted  with  some  success [47]. 

The  alternative  to  Reynolds  averaging  is  filtering.  The  filtering  process  is  to  spatially  filter  out  the 
turbulent  eddies  whose  scales  are  smaller  than  the  filter  width,  which  is  usually  mesh  size.  This  type  of 
approach  is  known  as  a  LES.  This  filtering  process  also  yields  new  terms  which  are  known  as  SGS  and 
must  be  modeled  to  provide  closure  to  the  set  of  governing  equations.  Namely,  large  eddies  are  resolved 
and  the  sub-grid  scale  eddies  are  modeled.  The  simplest  and  the  most  commonly  used  model  for  the  SGS 
is  proposed  by  Smagorinsky[48  ]  based  on  gradient-diffusion  concept  for  meteorological  application.  The 
eddy  viscosity,  which  represents  the  effects  of  turbulence,  is  calculated  from  an  algebraic  expression 
which  includes  a  model  constant,  the  modulus  of  the  rate  of  strain  tensor,  and  an  filter  width. 

Another  approach  is  DES  proposed  by  Spalart[49].  DES  employs  the  RANS  models  near  to  the  wall 
and  LES  in  the  wake  region  of  a  flow  where  unsteady  and  chaotic  motion  of  flow  is  usually  found. 
Namely,  RANS  model  is  used  away  from  the  wake  region  of  the  flow  to  save  computational  time 
compared  to  the  usage  of  LES,  and  LES  is  used  to  compute  the  eddies  and  vortices  to  maintain  the 
dynamic  features  of  the  flow  in  the  wake  region.  DES  has  been  studied  by  various  research  groups 
including  Boeing  Commercial  Airplane  group  and  ANSYS  CFX  group.  Boeing  Commercial  Airplane 
group  employs  one-equation  turbulence  model,  which  is  Spalart-Allmaras  model [49], while  ANSYS  CFX 
group  employs  Shear  Stress  Transport  (SST)  model[50]  in  the  RNAS  models  in  DES.  Here,  SST  model 
combines  the  positive  features  of  k  —  a  and  k  —  co  models.  In  other  words,  SST  employs  the  k  —  co 
model  near  wall  and  k  —  a  model  near  the  boundary  layer  edge. 
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II-3.  Solver  Selection 


Numerical  accuracy  and  computational  efficiency  are  the  key  criteria  to  select  the  flow  solver. 
Unsteady  flow  simulation  is  performed  at  the  Reynolds  number  of  1.0  X  104  for  the  flow  past  a  sphere 
to  decide  the  flow  solver.  Table  1  summarizes  Strouhal  numbers  and  drag  force  coefficients  (positive 
x  —axis  direction)  for  the  flow  past  a  sphere  at  a  Reynolds  number  of  1.0  X  104.  Strouhal  number  and 
drag  force  coefficient  are  defined  as  St  =  fx(D/V00 )  and  Cd  =  FD/(0.5pV^^4)respectively.  Here, 
f,  D,  Fd  Vqo,  p,  and  A  are  frequency  of  oscillating  drag  force  coefficient,  diameter  of  a  sphere,  drag 
force,  free-stream  velocity,  air  density,  and  cross-sectional  area  of  a  sphere,  respectively. 


Table  1.  Strouhal  number  and  drag  force  coefficient  at  Re  =  1.0  x  104 


Researcher  [Ref.] 

Tools 

Strouhal  Number 

Mean  CD 

Constantinescu 
and  Squires  [25] 

Calculation  (LES  model) 

0.20 

0.393 

Calculation  (DES  model) 

0.20 

0.397 

Jones  and  Clarke  [17] 

Calculation  (LES  model) 

0.19 

0.387 

Constantinescu,  ,Chapelet, 
and  Squires  [24] 

Calculation  (LES  model) 

0.20 

0.393 

Calculation  (DES  model) 

0.20 

0.397 

Achenbach  [27] 

Experiment 

N/A 

0.40 

Present 

Fluent  (LES  model) 

0.20 

0.451 

CFX  (LES  model) 

0.17 

0.445 

Table  1  compares  the  Strouhal  numbers  and  time  averaged  drag  force  coefficients.  In  general,  the 
Strouhal  numbers  for  the  flow  past  a  sphere  exist  between  0.19  and  0.20.  Fluent  gives  0.2,  on  the  other 
hand  CFX  does  0.17  as  shown  in  Fig.  II— 1.  In  Fig.  II— 1,  the  amplitude  and  Strouhal  number  of 
oscillating  drag  force  coefficient  are  obtained  through  the  FFT  (Fast  Fourier  Transform)  in  Fluent. 
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(a)  Fluent  (b)  CFX 

Fig.  II-l  Power  spectrum  of  the  drag  force  coefficient  at  fie  =  1.0  x  104 

Generally,  the  drag  force  coefficient  lies  between  0.39  and  0.40  for  the  flow  past  a  sphere.  Both 
Fluent  and  CFX,  however,  give  larger  drag  force  coefficient.  Figure  II — 2  shows  the  time  history  of  the 

TXV 

drag  force  coefficient.  Here  7*  represents  a  non-dimensional  time  which  is  defined  as  7*  =  d°°  and 

7  is  a  real  time  in  second. 


Fig.  II— 2  Time  history  of  the  drag  force  coefficient  at  Re  =  1. 0  x  104 

Table  2  compares  the  computational  speed.  Two  Intel  Quad-Core  Xeon  E5450,  3.0  GHz  CPU  are 
used  to  simulate  the  unsteady  flowfield  around  a  sphere  with  the  grid  size  of  2.2  million  cells.  It  seems 
that  Fluent  is  more  efficient  than  CFX  from  the  capability  of  parallel  computation  point  of  view  in  our 
simulations. 
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Table  2.  Comparisons  of  computational  speed 


Flow  solver 

Number  of  time  steps  per  day 

Fluent 

890  steps  (22,250  iterations) 

CFX 

266  steps  (3,993  iterations) 

Fluent  shows  reasonable  numerical  results  and  is  well  adapted  in  our  hardware  resources  so  that 
Fluent  is  selected  to  simulate  flowfield.  Here,  we  have  to  mention  that  even  though  we  select  Fluent  as 
our  flow  solver  for  the  present  research,  it  does  not  mean  that  Fluent  is  better  flow  solver  than  CFX 
officially.  Other  people  may  select  CFX  as  their  flow  solver. 

II-4.  Turbulence  Model  Selection 

There  are  many  kinds  of  different  turbulence  models  including  zero  equation  model,  one-equation 
model,  two-equation  models,  which  are  related  with  RANS,  LES,  and  DES. 

Jones  and  Clarke[17]  states  that  numerical  simulations  of  the  flow  past  a  sphere  have  been 
performed  mainly  using  LES  or  DES  in  the  range  of  Reynolds  numbers  greater  than  1.0  X  104.  They 
simulated  unsteady  flow  around  a  sphere  using  LES  at  the  Reynolds  numbers  of  1.0  X  104  and 
1.0  X  106.  Constantinescu  and  Squires [25]  also  says  that  DES  and  LES  showed  better  agreement  with 
measurements  than  URANS. 

Therefore,  LES  and  DES  are  considered  as  candidate  turbulence  models  in  the  present  research. 
Numerical  simulations  using  LES  and  DES  are  performed  for  the  flowfield  around  a  sphere  at  the 
Reynolds  number  of  1.0  X  104. 


(a)  LES  (b)  DES-SA 

Fig.  II— 3  Grid  structures  near  the  surface  (y+  =  1) 
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Figure  II— 3  shows  the  grid  structures  near  the  solid  wall,  which  are  used  for  the  current  LES  and 
DES-SA(Spalart-Allmaras)  simulations.  The  grid  density  near  the  solid  wall  for  a  given  Reynolds  number 
should  be  fine  enough  to  resolve  the  smaller  scale  eddies  that  arise  via  shear  layer  instabilities.  Here,  the 
grid  Reynolds  number  of  y+  =  1  is  used  for  the  very  first  grid  point  from  the  wall  to  have  the  best 
prediction  of  the  flowfield.  Here  the  grid  Reynolds  number  is  defined  as  y+  =  ^  and  y,  uT,  and  v 


represents  the  distance  from  the  wall  to  the  very  first  grid  point,  friction  velocity,  and  kinematic  viscosity 


respectively.  In  addition,  friction  velocity  is  defined  as  u. 


_  [l  =  v  Iel 
KoCW  2 


where  r  is  wall  shear  stress. 


When  the  skin  friction  coefficient  Cf  is  taken  as  ~  0.03 7ReD^S  from  the  empirical  correlations,  y 

_ 9_ 

becomes  y  =  y+  -j===ReD10 .  Grid  for  LES  is  prepared  using  Gambit  and  DES-SA  grid  is  generated 
based  on  the  Ref.  [25]. 


Figure  II— 4  shows  time  history  of  drag  force  coefficient.  DES  hires  Spalart-Allmaras  model[49]  in 
RANS  equations  and  seems  to  have  a  lot  of  noise  signal  compared  with  LES. 

Figure  II— 5  shows  Strouhal  number  of  drag  force  oscillations.  LES  shows  better  acceptable 
numerical  result  than  DES. 
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(a)  LES  (b)  DES-SA 

Fig.  II— 5  Power  spectrum  of  the  drag  force  coefficient  at  Re  =  1. 0  x  104 


(a)  LES  (b)  DES-SA 

Fig.  II-6  Visualization  of  vortex  structures  using  Q-criterion  at  Re  —  1. 0  x  104 


Figure  II— 6  plots  vortex  structures  near  the  sphere  using  the  second  invariant  of  the  velocity 
gradient [51]  called  Q-criterion  on  a  plane  out  of  the  three-dimensional  domain.  LES  shows  small  size 
eddies  near  the  surface  which  is  related  with  shear  layer  instability.  However  DES  with  Spalart-Allmaras 
turbulence  model  does  not  show  small  size  eddies  like  LES.  Actually,  the  flowfield  around  a  sphere  at 
Re  =  104  has  turbulent  boundary  layer  on  the  surface  and  turbulent  wake  region  in  the  base  region[23]. 
LES,  therefore,  gives  better  numerical  results  than  DES  and  is  selected  as  a  turbulence  model  in  the 
present  research. 

II-5.  Solver  Validation 

Numerical  validation  for  Fluent  has  been  performed  for  the  flow  past  a  sphere  at  the  Reynolds 
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numbers  of  1.0  x  103,  2.0  x  103,  4.0  x  103,  and  1.0  x  104.  Sakamoto  and  Haniu[23]  states  that  the 
flow  structure  behind  a  sphere  is  turbulent  wake  and  the  flow  along  the  surface  is  laminar  at  the  Reynolds 
numbers  between  800  and  3.0  X  103.  When  the  Reynolds  number  exceeds  3.0  X  103,  the  vortex  sheet 
which  separates  from  the  sphere  surface  changes  from  laminar  to  turbulent  flow.  In  addition,  the  vortex 
sheet  separating  from  the  surface  of  the  sphere  becomes  completely  turbulent  when  the  Reynolds  number 
exceeds  6.0  X  103.  At  the  Reynolds  number  ranging  from  8.0  X  102  to  1.5  X  104,  higher  and  lower 
frequency  modes  are  coexisted  because  periodic  fluctuation  in  the  vortex  tube  formed  by  the  pulsation  of 
the  vortex  sheet  separated  from  the  surface  of  the  sphere  and  in  the  turbulent  waving  wake. 

Steady-state  flows  past  a  sphere  of  diameter  D  =  0.1  m  are  obtained  by  solving  RANS  equations 
with  Spalart-Allmaras  turbulence  model  in  the  following  computational  domain  whose  far-field  size  is  50 
times  of  sphere  diameter  as  shown  in  Fig.  II— 7.  The  density  and  viscosity  of  the  free-stream  is  p  = 
1.225 kg/m3  and  /r  =  1.7894  x  10 ~skg/(m  •  s)  ,  respectively.  The  free-stream  velocity  can  be 
calculated  by  Reynolds  number,  for  example  Vk>  =  1.46073m/s  at  Re  =  1.0  x  104.  The  computational 
domain  has  the  velocity-inlet  with  0.1%  of  turbulent  intensity,  pressure-outlet,  and  no-slip  boundary 
condition  on  the  surface  of  a  sphere.  The  number  of  grid  cells  is  2,212,656  to  2,263,845.  Two  CPUs  of 
Intel  Zeon  X5450  3.0  GHz  are  used  to  compute  the  steady-state  flow  solutions. 


(a)  Shape  of  computational  domain 


(b)  Far-field  size  of  computational  domain 


Fig.  II-7  Computational  domain  of  flowfield  past  a  sphere 


Unsteady  flowfield  for  LES  simulation  is  initialized  using  steady-state  flow  condition.  A  physical 
time  step  is  chosen  as  0.02D/I4,  [24],  which  insures  approximately  250  time  steps  per  cycle, 
corresponding  to  the  main  shedding  frequency,  for  example,  in  case  of  Re  =  1.0  X  104.  Each  physical 
time  step  has  25  sub-iterations  to  have  fully  converged  solution.  Two  CPUs  of  Intel  Zeon  X5450  3.0  GHz 
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spend  about  2  to  2.3  minutes  per  time  step  and  each  LES  simulation  has  10,000  to  13,000  time  steps. 
Usually,  15  to  20  days  are  required  to  finish  one  unsteady  flow  simulation. 


Fig.  II-8  Time  histories  of  drag  force  coefficients  at  different  Reynolds  numbers 

Figure  II — 8  compares  time  history  of  drag  force  coefficients  at  different  Reynolds  numbers  of 
1.0  X  103,  2.0  X  103,  4.0  X  103,  and  1.0  X  104.  Table  3  summarizes  drag  force  coefficients  again. 
Here,  mean  drag  force  coefficient  for  unsteady  flow  is  obtained  by  taking  time  average  of  unsteady  drag 
force  coefficients  and  drag  force  coefficient  for  steady-state  flow  is  obtained  by  solving  steady-state  flows. 
Unsteady  flows  using  LES  predict  less  drag  force  coefficient  than  steady-state  flows  using  RANS  with 
Spalart-Allmaras  turbulence  model.  Figure  II— 9  plots  the  drag  force  coefficient  along  a  smooth  sphere  as 
a  function  of  Reynolds  number  and  confirms  that  the  present  numerical  results  are  very  reasonable.  Jones 
and  Clarke[17]  also  states  that  the  drag  force  coefficients  are  close  to  0.5  in  the  Reynolds  number  range  of 
800  to  3.8  x  105. 


Table  3.  Mean  drag  force  coefficients  for  a  sphere  at  different  Reynolds  numbers 


Reynolds  Number 

Mean  Cd  (Unsteady) 

Cd  (Steady-state) 

1,000 

0.513 

0.515 

2,000 

0.468 

0.486 

4.000 

0.444 

0.470 

10,000 

0.451 

0.454 

- 19- 


Amplitude  Amplitude 


Fig.  II-9  Drag  force  coefficient  of  a  smooth  sphere  as  a  function  of  Reynolds  number 


(a)  Re  =  1.0  x  103  (b)  Re  =  2.0  x  103 


(c)  Re  =  4.0  x  103  (d)  Re  =  1.0  x  104 

Fig.  11-10  Power  spectrum  of  the  drag  force  coefficient 
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Strouhal  Number 


(a)  Present  numerical  data 


(b)  Available  numerical/experimental  data  (Ref.[23],  pp.388,  Fig.  3) 

Fig.  11-11  Various  data  for  Strouhal  number  of  sphere 
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Figure  II- 10  plots  the  power  spectrum  of  the  drag  force  coefficient  in  association  with  Fig.  II— 8  at 
different  Reynolds  numbers.  High  frequency  mode  due  to  vortex  sheet  separating  from  the  surface,  which 
is  called  Kelvin-Helmholtz  instability,  as  well  as  low  frequency  mode  due  to  large  scale  vortex  shedding 
is  observed  although  high  frequency  mode  is  not  detected  at  the  Reynolds  number  of  1.0  X  104  in  the 
present  simulation.  Sakamoto  and  Haniu[23]  also  states  that  Cometta  obtained  high  frequency  mode  only 
up  to  Re  =  7.4xl03  in  their  paper.  Even  though  the  Strouhal  numbers  differ  from  researcher  to 
researcher  in  the  range  of  103  <  Re  <  104[23],  the  present  Strouhal  numbers  look  very  reasonable 
compared  to  Fig.  II—  1 1 . 


(a)  flow  visualization  past  a  sphere  using  particle  trace 


(b)  patterns  of  vortex  shedding  in  the  wake  region  using  Q-criterion 

Fig.  11-12  Flow  visualization  past  a  sphere  at  Re  —  2. 0  x  103 

Figure  11-12  visualizes  the  flow  patterns  past  a  sphere  at  the  Reynolds  number  of  2.0  X  103  using 
particle  trace  technique  and  Q-criterion.  The  pattern  of  vortex  shedding  in  the  wake  region  is  very  similar 
to  the  pattern  of  the  particle  trace  in  that  region  so  that  the  flow  visualization  technique  using  Q-criterion 
can  replace  the  flow  visualization  using  particle  trace. 

Figure  11-13  plots  mean  skin  friction  coefficient,  Cf ,  distributions  along  the  surface  of  a  sphere  in 
the  streamwise  direction.  There  is  a  reverse  flow  in  the  streamwise  direction  and  flow  separation  occurs 
where  the  positive  skin  friction  coefficient  changes  to  the  negative  skin  friction  coefficient  or  vice  versa. 
Table  4  tabulates  the  separation  locations  according  to  the  Fig.  11-13.  Lower  polar  angles  represent  the 
locations  of  separation  on  the  top  surface  and  higher  ones  represent  the  locations  of  separation  on  the 
bottom  surface. 

Jones  and  Clarke[17]  estimated  that  separation  occurs  at  a  polar  angle  of  88.0°  ±  1°  at  the 
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Reynolds  number  of  1.0  X  104.  At  the  same  Reynolds  number,  Constantinescu  et  al.[24]  obtained 
85.0°  ±  1°  using  LES.  The  present  simulation  using  LES  obtained  86.4°  —  88.2°  on  the  top  surface  and 
270.0°  —  271.8°  on  the  bottom  surface  at  the  Reynolds  number  of  1.0  X  104  as  summarized  in  Table  4 
and  has  a  good  agreement  with  other  numerical  results. 


Fig.  11-13  Mean  skin  friction  coefficient  distributions  along  the  streamwise  direction 


Table  4.  Location  of  flow  separation  line  as  polar  angles 


Reynolds  Number 

Flow  Separation  Angle  (degree) 

1,000 

99.0-100.8 

257.4-259.2 

2,000 

93.6-95.4 

264.6-266.4 

4,000 

90.0-91.8 

268.2-270.0 

10,000 

86.4-88.2 

270.0-271.8 

Figure  11-14  plots  the  oil  flows  on  the  top  and  bottom  surface  of  the  sphere  at  different  Reynolds 
numbers  with  the  different  legend  of  skin  friction  coefficient. 

As  mentioned  in  Ref.  [23],  the  vortex  sheet  which  separates  from  the  sphere  surface  is  laminar  flow 
as  shown  in  Fig.  II-14(a)&(b).  At  the  Reynolds  number  of  4.0  X  103,  the  vortex  sheet  changes  from 
laminar  to  turbulent  flow  as  shown  in  Fig.  II-14(c).  Finally,  the  vortex  sheet  becomes  completely 
turbulent  flow  as  shown  in  Fig.  II— 14(d)  so  that  the  band  shape  of  flood  type  contours  for  skin  friction 
coefficient  becomes  wiggled  along  the  latitude  of  a  sphere.  The  flow  separation  lines  are  also  observed 
along  the  surface. 
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(a)  Re  =  1.0  x  103 


(b)  Re  =  2.0  x  103 


(c)  Re  =  4.0  x  103 


(d)  Re  =  1.0  x  104 

Fig.  11-14  Oil  flow  on  the  top  and  bottom  surface  at  certain  instantaneous  time 
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III.  NUMERICAL  RESULTS  AND  DISCUSSIONS 


As  mentioned  in  the  statement  of  work,  the  objective  of  this  research  is  to  figure  out  the  effects  of 
geometric  shapes  and  Reynolds  numbers  on  the  unsteady  flow  physics  including  vortex  shedding 
formation,  drag,  lift,  and  side  forces  exerted  on  the  body  to  prepare  a  design  guideline  for  MAV. 
Rectangular  parallelepipeds  and  circular  cylinders  are  therefore,  prepared  as  shown  in  Fig.  Ill— 1.  Each 
geometry  has  three  different  aspect  ratios  of  1.0,  1.5,  and  2.0  as  shown  in  Fig.  Ill- 1 .  Sphere  has  always 
one  aspect  ratio  of  1.0.  Analysis  of  unsteady  flowfield  around  these  blunt  bodies  including  a  sphere  have 
been  accomplished  at  the  Reynolds  numbers  of  1.0  X  103,  2.0  X  103,  4.0  X  103  and  1.0  X  104. 


Rectangular  parallelepiped  (D  =  H  =  0.1  m) 

W/H=1.0 

W/H=1.5 

W/H=2.0 

H 

t 

H 

H 

'  W'  D 

Circular  cylinder  (D  =  0.1  m) 

W/D=1.0 

W/D=1.5 

W/D=2.0 

D 

D 

V 

Fig.  Ill— 1  Blunt  body  geometries  for  unsteady  flow  simulation 

The  present  simulations  are  performed  in  the  rectangular  parallelepiped  type  computational  domain 
whose  size  is  40D  X  40 D  X  75 D  as  shown  in  Fig.  Ill — 2,  which  has  the  velocity-inlet,  pressure-outlet, 
and  symmetry  conditions  on  the  top,  bottom,  and  two  lateral  sides.  In  Figs.  III-2(b)&(c),  the  right  side 
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view  means  the  view  from  the  positive  z -direction  to  the  negative  z -direction.  No-slip  boundary 
condition  is  imposed  on  the  surface  of  the  body.  On  the  other  hand,  the  computational  domain  and 
boundary  conditions  shown  in  Fig.  II— 7  are  applied  to  the  simulation  of  the  flowfield  past  a  sphere. 

Steady-state  flow  condition  is  used  as  an  initial  flow  condition  for  unsteady  flow  simulations.  The 
density  and  viscosity  of  the  free-stream  are  p  =  1.22Skg/m3  and  fi  —  1.7894  x  10~skg/(m  •  s)  , 
respectively.  The  free-stream  velocity  can  be  calculated  by  Reynolds  number  definition.  A  physical  time 
step  is  chosen  as  0.02D/14>  and  each  time  step  has  25  sub-iterations  to  make  converged  status. 


(c)  Right  side  view  of  circular  cylinder  (c)  Top  view  of  rectangular  parallelepiped/circular  cylinder 

Fig.  III-2  Computational  domain 

III-l.  Analysis  of  Flowfield  Past  a  Sphere 
1.  Vortex  shedding  formation 

Figure  III— 3  shows  the  vortex  distributions  around  a  sphere  using  Q-criterion  at  the  Reynolds 
number  of  1.0  X  104  after  the  vortex  sheds  from  the  top  surface. 
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(a)  Right  side  view 


(b)  Top  view 

Fig.  III-3  Visualization  of  a  hairpin-shaped  vortex  shedding  around  a  sphere  at  Re  =  1. 0  x  104 


(b)  Velocity  distributions  (Right  side  view) 


(d)  Velocity  distributions  (Top  view) 


Fig.  III-4  Flow  features  at  Re  =  1. 0  x  104 
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Figure  III— 4  plots  the  detailed  flow  features  near  a  sphere  at  the  Reynolds  number  of  1.0  X  104  to 
understand  the  vortex  formation  process.  The  detached  vortex  loop  from  the  surface  flows  to  the 
downstream  direction.  The  upper  part  of  the  vortex  loop  flows  faster  than  the  lower  part  of  the  vortex  loop 
and  spreads  out  while  making  a  hairpin-shaped  vortex  tube.  The  lower  part  of  the  vortex  loop  flows 
slowly  to  the  downstream  direction  while  making  two  long  legs  because  of  recirculating  flow  to  the  body 
as  shown  in  Fig.  Ill — 4(b)  and  inward  directional  flow  to  the  center  of  the  body  as  shown  in  Fig.  Ill-  4(d). 


Fig.  III-5  Flow  visualization  near  a  sphere  at  Re  =  1. 0  x  103 

Figure  III— 5  repeats  the  above  description  again  at  a  different  Reynolds  number.  The  vortex  loop 
detached  from  the  body  flows  to  the  downstream  direction  with  different  flow  speed  conditions;  upper 
part  of  the  vortex  loop  flows  faster  and  lower  part  of  the  vortex  loop  flows  slower,  so  that  the  vortex  loop 
leans  as  shown  in  Fig.  Ill— 5 .  The  structure  of  the  lower  part  of  the  vortex  loop  is  very  well  demonstrated 
by  the  stream  lines. 
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(a)  Velocity  distributions  of  u-component  (right  side  view) 
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(b)  Velocity  distributions  of  w-component  (bottom  view) 

Fig.  III-6  Velocity  component  distributions  on  the  vortex  structure  at  Re  =  1. 0  x  104 


The  streamwise  directional  velocity  components  (u  -component)  and  the  transverse  directional 
velocity  components  (w-component)  are  distributed  on  the  vortex  structures  as  shown  in  Fig.  III-6.  As 
explained  above,  the  upper  part  of  the  vortex  loop  (hairpin-shaped  vortex  tube  colored  by  orange)  flows 
to  the  downstream  direction  faster  than  the  lower  part  of  the  vortex  loop  (two  long  legs  colored  by  green) 
as  shown  in  Fig.  III-6(a).  The  distance  between  two  long  legs  becomes  narrow  while  fluid  flows  to  the 
downstream  direction  because  there  are  inward  directional  flows  (green  and  sky  blue)  to  the  center  of  the 
sphere  as  shown  in  Fig.  III-6(b). 

The  vortex  shedding  occurs  from  one  position  to  another  position  along  the  circumference  of  a 
sphere  and  results  in  the  oscillating  drag,  lift,  and  side  forces  as  shown  in  Fig.  Ill— 7.  Here,  Cd ,  Q,  and 
Cs  represent  the  drag,  lift,  and  side  force  coefficient,  respectively. 


Time  history  of  unsteady  force  coefficient  at/?e  =  1.0xl04 
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Fig.  III-7 


The  location  of  the  vortex  shedding  changes  from  time  to  time.  Figure  III— 8  illustrates  the  irregular 
rotation  of  the  vortex  shedding  location  about  streamwise  axis  through  the  center  of  the  sphere.  Pressure 
contours  are  plotted  on  the  y  —  z  plane  at  different  time  status  shown  in  Fig.  III-7,  where  it  locates  at 
x/D  =  0.6  from  the  center  of  the  sphere.  The  lower  and  higher  pressure  region  is  rotating  as  time  elapses. 


P- P' P' P- £■  $>■  £■  P- P- 


Pressure  [Pa] 

Fig.  III-8  Pressure  contours  at  the  location  of  x/D  —  0. 6  at  Re  —  1.0  x  104 
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Figure  III— 9  shows  vortex  sheets  separating  from  the  surface  of  a  sphere  at  the  Reynolds  number  of 
1.0  X  103.  The  high  mode  Strouhal  number  of  oscillating  force  coefficient  comes  from  the  vortex  sheet 
separation  from  the  surface. 

Taneda[22]  sketched  vortex  sheet  separating  from  a  sphere  in  the  range  of  Reynolds  numbers 
between  3.8  X  105  and  1.0  X  106  as  shown  in  Fig.  III-10(a).  He  described  that  “the  vortex  sheet  shed 
from  the  sphere  rolls  up  to  form  a  pair  of  streamwise  vortices.”  Figure  III— 10(b)  plots  the  vortex  sheet 
separating  from  the  surface  using  Q-criterion  and  agrees  well  with  the  Taneda’s  sketch. 


(a)  sketch  of  the  vortex  sheet  [ref.  22]  (b)  vortex  sheet  using  Q-criterion 

Fig.  III-10  Vortex  sheet  separating  from  the  surface 


2.  Parametric  study 


Fig.  III-ll  Strouhal  number  versus  Reynolds  number 
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Figure  III- 11  plots  the  Strouhal  number  versus  Reynolds  number  for  the  sphere.  The  present 
numerical  result  shows  that  the  Strouhal  number  of  drag  force  coefficient  is  0.2  more  or  less  over  the 
entire  range  of  Reynolds  numbers  and  agrees  very  well  with  the  available  numerical  or  experimental 
results[17,  23].  The  high  mode  Strouhal  number,  which  is  associated  with  the  small  scale  instability  of  the 
separating  shear  layer  on  the  surface  is  also  very  reasonable.  However,  the  high  mode  Strouhal  number  at 
Re  =  1.0  x  104  was  not  able  to  obtain  in  the  present  simulation. 

The  Strouhal  number  of  the  lift  or  side  force  coefficient  is  less  than  that  of  drag  force  coefficient. 
Actually,  the  Strouhal  number  of  the  lift  force  is  half  of  Strouhal  number  of  the  drag  force  in  a  two- 
dimension  circular  cylinder  problem.  The  existence  of  Strouhal  number  for  the  lift  and  side  force 
coefficient  means  that  the  vortex  loops  continuously  shed  into  the  streamwise  direction  in  an  irregularly 
rotating  manner  with  respect  to  the  axis  parallel  to  the  downstream  through  the  center  of  the  sphere. 

The  high  mode  Strouhal  number  exists  only  in  the  streamwise  direction  and  strongly  depends  on  the 
Reynolds  number.  All  the  low  mode  Strouhal  numbers,  however,  almost  do  not  dependent  on  the 
Reynolds  numbers. 


Figure  III- 12  plots  the  amplitude  variations  according  to  the  variations  of  Reynolds  numbers.  The 
amplitude  of  drag  force  coefficient  is  very  small  compared  with  that  of  the  lift  or  side  force  coefficient.  In 
the  case  of  vortex  shedding  from  a  sphere,  Figs.  Ill— 3  and  III— 4  state  that  the  vortex  loop  detaches  from  a 
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sphere  rotates  irregularly  at  the  Reynolds  numbers  of  1.0  x  103,  2.  0  x  103,  4.0  x  103,  and  1.0  x  104. 
The  amplitude  of  lift  and  side  force  coefficient  increases  as  Reynolds  number  increases  although  the 
oscillating  frequency  is  almost  same.  The  amplitude  of  the  drag  force  coefficient  of  high  frequency  mode 
is  nearly  zero  even  though  the  oscillating  frequency  is  very  high. 

In  general,  the  amplitude  of  the  lift  force  or  side  force  coefficient  of  the  low  frequency  mode  is  much 
greater  than  that  of  drag  force  coefficient  of  the  low  frequency  mode  as  shown  in  Fig.  Ill— 12.  In  addition, 
the  amplitude  of  the  drag  force  coefficient  of  low  frequency  mode  is  greater  than  that  of  the  drag  force 
coefficient  of  high  frequency  mode  .  If  there  exists  a  high  frequency  mode  of  lift  or  side  force  coefficient, 
the  amplitude  of  the  lift  force  or  side  force  coefficient  of  high  frequency  mode,  which  is  directly  related 
with  the  amplitude  of  the  drag  force  coefficient  of  high  frequency  mode,  is  much  smaller  than  that  of  the 
lift  force  or  side  force  coefficient  of  low  frequency  mode  so  that  the  high  mode  Strouhal  number  of  the  lift 
force  or  side  force  coefficient  cannot  be  observed  in  FFT.  In  other  words,  the  high  mode  Strouhal  number 
of  the  lift  force  or  side  force  coefficient  can  not  be  plotted  in  Fig.  Ill— 11. 


Fig.  III-13  Mean  force  coefficient  versus  Reynolds  number 


Figure  III- 13  plots  the  mean  force  coefficient  versus  Reynolds  number.  As  already  mentioned  in 
section  II-5,  the  mean  drag  force  coefficient  is  about  0.5  at  the  Reynolds  numbers  from  1.0  X  103to 
1.0  X  104.  The  vortex  shedding  from  a  sphere  does  not  occur  in  an  axisymmetric  pattern.  Actually,  the 
vortex  sheds  form  the  body  in  randomly  rotating  pattern  about  an  axis  parallel  to  the  downstream  through 
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the  center  of  a  sphere  so  that  mean  lift  force  coefficient  and  side  force  coefficient  for  the  flow  past  a 
sphere  may  not  be  zero  as  shown  in  Fig.  Ill—  1 3  even  though  they  are  almost  zero  values. 


z/D 

3- 


x=3.2D 


0.0  0.5  1.0 

m/VL 


z/D 

3- 


x=6.4D 


0.0  0.5  1.0 

u/V„ 


Fig.  III-14  Velocity  profiles  of  time  averaged  streamwise  velocity  component 


Figure  III— 14  plots  the  time  averaged  streamwise  velocity  profiles  at  several  different  distances  from 
the  body.  The  time  averaged  streamwise  velocity  component  is  normalized  based  on  the  free-stream 
velocity  and  plotted  along  the  z-direction  at  the  middle  section  of  the  height.  The  recirculation  region 
disappears  faster  as  the  Reynolds  number  increases. 


3.  Summary 

The  formation  mechanism  of  the  hairpin-shaped  vortex  shedding  and  the  unsteady  oscillating  flow 
phenomena  are  numerically  investigated. 

There  are  two  mode  Strouhal  numbers  in  the  range  of  Reynolds  numbers  from  1.0  X  103  to 
1.0  X  104;  one  is  the  low  mode  Strouhal  number  which  comes  from  the  large  scale  vortex  shedding  and 
the  other  is  high  mode  Strouhal  number  which  comes  from  the  Kelvin-Helmholtz  instability  due  to  vortex 
sheet  separating  from  the  body  surface.  The  high  mode  Strouhal  number  exists  only  in  the  streamwise 
direction  and  strongly  depends  on  the  Reynolds  number.  The  low  mode  Strouhal  numbers,  however, 
almost  do  not  depend  on  the  Reynolds  numbers.  The  vortex  sheds  from  the  body  in  randomly  rotating 
pattern  about  an  axis  parallel  to  the  downstream  direction  through  the  center  of  a  sphere. 

The  Strouhal  number  of  the  drag  force  coefficient  is  greater  than  that  of  lift  or  side  force  coefficient. 
The  amplitude  of  the  drag  force  coefficient  is  very  small  compared  with  that  of  the  lift  or  side  force 
coefficient.  The  amplitude  of  the  drag  force  coefficient  of  high  frequency  mode  is  nearly  zero  even 
though  the  oscillating  frequency  is  very  high. 

The  amplitude  of  lift  and  side  force  coefficient  increases  as  Reynolds  number  increases  although  the 
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oscillating  frequency  is  almost  same. 

Finally,  we  can  conclude  that  there  are  two  frequency  modes  of  Strouhal  numbers  for  the  drag  force 
coefficient;  one  is  the  low  mode  Strouhal  number  and  the  other  is  the  high  mode  Strouhal  number.  In 
addition,  the  location  of  the  vortex  shedding  rotates  irregularly  about  an  axis  through  the  center  of  the 
sphere.  The  Reynolds  number  strongly  affects  the  high  mode  Strouhal  number  of  the  drag  force 
coefficient  and  the  amplitudes  of  the  lift  and  side  force  coefficient. 

III-2.  Analysis  of  Flowfield  Past  a  Rectangular  Parallelepiped 
1.  Vortex  shedding  formation 

The  formation  process  of  the  hairpin-shaped  vortex  loop  is  investigated  using  Q-criterion  at  the 
Reynolds  number  of  1.0  X  103.  Figure  III- 15  shows  the  feature  of  vortex  shedding  from  a  cube  whose 
aspect  ratio  is  equal  to  1.0.  Figure  III- 16  shows  the  instantaneous  screen  shot  right  after  a  top  vortex 
shedding  and  at  the  very  beginning  stage  of  bottom  vortex  shedding. 


Fig.  III-15  Visualization  of  hairpin-shaped  vortex  shedding  using  Q-criterion 

(A/R=1.0,  Re  =  1.0  x  103) 


Fig.  III-16  Right  side  view  of  vortex  shedding  near  the  body  (A/R= 1.0,  Re  —  l.Ox  103) 

The  stage  in  Fig.  Ill—  1 6  can  be  explained  in  detail  as  follows:  The  separated  flows  at  the  front  edge 


on  the  top  surface  flow  over  the  surface  to  the  downstream  direction  and  part  of  them  recirculate  toward 
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the  front  edge  along  the  surface  then  they  flow  to  the  downstream  direction  with  the  separated  flows  at  the 
front  edge  together  and  are  detached  from  the  body.  Figure  III- 17  shows  the  detailed  flow  patterns  near 
the  surface  when  vortex  sheds  from  the  top  surface. 


(a)  Recirculation  flows  along  the  surface 


(b)  Vortex  flows  on  and  near  the  surface 


(c)  Vortex  distributions  (front  view) 


(d)  Vortex  distributions  (rear  view) 


(e)  Streamlines  near  the  top  surface 

Fig.  III-17  Vortex  shedding  on  the  top  surface  ( A/R=1.0 ,  Re  =  1.  0  x  103) 

After  the  vortex  loop  sheds  from  the  body,  the  high-velocity  side  of  vortex  shedding  (upper  part) 
flows  faster  than  the  low-velocity  side  of  vortex  shedding  (lower  part)  and  spreads  out  while  making 
round  tube  shape  like  a  hairpin  as  shown  in  Fig.  Ill- 18.  The  low-velocity  side  of  detached  vortex  loop 
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flows  to  the  downstream  direction  slower  than  the  high-velocity  side  of  vortex  shedding  with  the  shape  of 
two  long  legs. 


(a)  Vortex  distribution  near  the  body 


(b)  Velocity  distributions  near  the  body 

Fig.  III-18  Flow  features  near  the  body  (right  side  view  at  middle  section) 

(A/R=1.0,  Re  —  1.0  x  103) 


Fig.  III-19  Top  view  of  hairpin-shaped  vortex  (A/R=1.0,  Re  —  1. 0  x  103) 
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Figure  III- 19  shows  the  vortex  shedding  patterns  in  the  top  view.  While  the  upper  part  of  the  vortex 
shedding  spreads  out,  the  distance  between  the  two  long  legs,  which  is  the  lower  part  of  the  vortex 
shedding,  becomes  narrow  because  of  inward  directional  flow  as  shown  in  Fig.  Ill— 20. 


(a)  Velocity  distributions  near  the  body 


(b)  Vortex  distributions  near  the  body 


(c)  Vortex  distributions  near  the  body  (transparency  mode) 

Fig.  III-20  Flow  features  near  the  body  (top  view  at  middle  section)  (A/R=1.0,  Re  =  1. 0  x  103) 
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(a)  Velocity  distributions  of  u-component  (left  side  view) 
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(b)  Velocity  distributions  of  w-component  (top  view) 

Fig.  III-21  Velocity  component  distributions  on  the  vortex  structure  (A/R= 1.0,  Re  =  1. 0  x  104) 


Figure  III-21  demonstrates  the  mechanism  of  vortex  shedding  formation  in  another  way  using  Q- 
criterion  after  the  vortex  sheds  from  the  bottom  surface.  Two  velocity  components  of  u  and  w  are 
distributed  on  the  vortex  structures.  The  hairpin-shaped  vortex  tube  (lower  part  of  the  vortex  loop  colored 
by  orange)  flows  faster  to  the  downstream  direction  whereas  the  two  long  legs-shaped  vortex  tube  (upper 
part  of  the  vortex  loop  colored  by  light  green)  flows  slowly  to  the  downstream  direction  as  shown  in  Fig. 
III-21(a).  Figure  III-21(b)  demonstrates  that  the  distance  between  the  two  long  legs  (colored  by  sky  blue 
and  yellow)  decreases  because  of  inward  directional  flows  to  the  center  of  the  body. 
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(a)  Velocity  distributions  of  u-component  (right  side  view) 

V**ocity  w 


(b)  Velocity  distributions  of  w-component  (bottom  view) 

Fig.  III-22  Velocity  component  distributions  on  the  vortex  structure  (A/R= 2.0,  Re  =  1. 0  x  104) 


Figure  III-22  shows  another  simulation  result.  Figure  III-22(b)  explains  that  the  hairpin-shaped 
vortex  tube  (left  circle  in  the  figure)  spreads  out  and  the  two  long  legs  (right  circle  in  the  figure)  approach 
to  each  other. 


Fig.  III-23  Time  history  of  unsteady  force  coefficient  ( A/R=1.0 ,  Re  =  1. 0  x  104) 
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(a)  plane  location 


(c)  T*  =  110 


(b)  T*  =  70 
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(d)  r  =  130 


Fig.  III-24 
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Pressure  [Pa] 

Pressure  contours  at  the  location  of  x/D  =  0. 6  G4//?  =1.0,  Re  =  1. 0  x  104) 


Figure  III-24  explains  the  rotating  tendency  using  the  pressure  contours  on  the  y  —  z  plane  at 
different  time  status  shown  in  Fig.  III-23  at  the  aspect  ratio  of  1.0.  The  pressure  contours  changes  about 
an  axis  through  the  center  of  the  body  from  time  to  time,  which  means  that  the  location  of  vortex 
shedding  rotates  about  an  axis.  Figures  III-25&III-26  explain  the  same  phenomena  at  the  aspect  ratio  of 
2.0. 
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Time  (T‘) 

Fig.  III-25  Time  history  of  unsteady  force  coefficient  (A/R-2A),  Re  —  1. 0  x  104) 
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Fig.  III-26  Pressure  contours  at  the  location  of  x  x/D  =  0. 6  (A/R  =2.0,  fie  =  1.0x  104) 
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Fig.  III-27  Sequence  of  hairpin-shaped  vortex  formation 
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Figure  III— 27  shows  the  sequence  of  hairpin-shaped  vortex  formation  in  certain  time  period  and 
states  that  the  hairpin-shaped  vortex  rotates  as  time  elapses. 

2.  Parametric  study 


Fig.  III-28  Strouhal  number  versus  Reynolds  number 

Figure  III— 28  plots  the  Strouhal  number  versus  Reynolds  number  for  the  rectangular  parallelepiped 
with  different  aspect  ratios.  The  variation  of  Strouhal  number  according  to  the  variation  of  Reynolds 
number  and  aspect  ratio  looks  quite  complicated.  As  might  be  expected,  the  Strouhal  number  of  drag 
force  coefficient  is  higher  than  that  of  lift  force  coefficient  or  side  force  coefficient  for  the  three  different 
aspect  ratios.  The  lift  and  side  force  coefficients  are  defined  as  the  same  manner  as  the  drag  force 
coefficient  whose  reference  area  is  D  XW.  In  addition,  the  Strouhal  number  of  drag  force  coefficient 
strongly  depends  on  the  Reynolds  number  as  well  as  the  aspect  ratio.  In  case  of  aspect  ratio  of  1.0,  the 
Strouhal  numbers  for  lift  and  side  force  coefficient  can  be  interchangeable  because  of  geometric  shape.  In 
other  words,  the  Strouhal  number  of  lift  force  coefficient  can  be  regarded  as  the  Strouhal  number  of  side 
force  coefficient  or  vice  versa.  However,  the  Strouhal  numbers  for  the  lift  and  side  force  coefficient 
cannot  be  interchangeable  when  the  aspect  ratio  is  1.5  or  2.0  because  the  front  section  of  the  geometry  is 
not  a  square  anymore.  The  Strouhal  number  of  the  lift  force  coefficient  is  higher  than  the  Strouhal  number 
of  the  side  force  coefficient  in  cases  of  the  aspect  ratios  of  1.5  and  2.0  because  the  vortex  shedding  from 
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edge  of  the  width  is  more  dominant  than  the  vortex  shedding  from  the  edge  of  the  height.  In  other  words, 
vortex  loop  detaches  more  frequently  from  the  edge  of  the  width  than  from  the  edge  of  the  height.  In  case 
of  aspect  ratio  of  1.0,  vortex  loop  almost  evenly  detaches  from  the  edge  of  the  width  and  edge  of  the 
height. 


Figure  III— 29  plots  the  amplitude  variations  of  the  force  coefficients  according  to  the  variations  of 
Reynolds  numbers.  The  amplitude  of  drag  force  coefficient  is  very  small  compared  with  that  of  the  lift  or 
side  force  coefficient.  The  amplitudes  of  both  lift  and  side  force  coefficients  increase  as  Reynolds  number 
increases  at  the  aspect  ratio  of  1.0.  The  amplitude  of  lift  force  coefficient  is  almost  same  magnitude 
except  at  the  Reynolds  number  of  2.0  X  103  meanwhile  the  amplitude  of  side  force  coefficient  decreases 
as  Reynolds  number  increases  when  the  aspect  ratio  is  1.5.  At  the  aspect  ratio  of  2.0,  the  amplitude  of  lift 
force  coefficient  increases  as  Reynolds  number  increases  while  the  amplitude  of  side  force  coefficient  is 
almost  same  magnitude  as  Reynolds  number  increases. 

In  the  range  of  Reynolds  numbers  1.0  X  103,  2.0  X  103,  and  4.0  X  103,  the  amplitude  of  lift  force 
coefficient  as  well  as  the  amplitude  of  side  force  coefficient  is  scattered  at  the  same  Reynolds  number  for 
the  different  aspect  ratios.  However,  the  amplitude  of  lift  force  coefficient  is  almost  same  magnitude  at 
the  Reynolds  number  of  1.0  X  104  no  matter  what  the  aspect  ratio  is  and  the  amplitude  of  side  force 
coefficient  is  also  almost  same  magnitude  at  this  Reynolds  number  except  the  case  of  aspect  ratio  of  1.0. 
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Figures  III— 28  and  III — 29  can  deduce  that  the  vortex  shedding  at  the  aspect  ratio  of  1.0  equally 
affects  the  lift  force  coefficient  oscillation  and  side  force  coefficient  oscillation,  and  sheds  in  an  irregular 
rotating  manner,  which  means  the  location  of  the  vortex  shedding  rotates  about  an  axis  through  the  center 
of  the  body.  In  cases  of  aspect  ratios  of  1.5  and  2.0,  the  rotating  tendency  becomes  weaker  as  the 
Reynolds  number  increases  as  shown  in  Fig.  Ill— 30  whereas  it  becomes  stronger  at  the  ratio  of  1.0.  In 
addition,  rotating  tendency  becomes  weaker  as  the  aspect  ratio  increases  as  shown  in  Fig.  Ill— 3 1 . 


(a)  Right  side  view  (Re  —  1.0  x  103) 


(c)  Right  side  view  (Re  —  1.0  x  104) 


(b)  Top  view  (Re  —  1.0  x  103)  (d)  Top  view  (Re  —  1.0  x  104) 

Fig.  III-30  Vortex  shedding  patterns  at  different  Reynolds  numbers  (A/R  —  1. 5) 


(a)  Right  side  view  (A/R  =  1.0) 


(c)  Right  side  view  (A/R  =  2.0) 


(b)  Top  view  (A/R  =  1.0) 


(d)  Top  view  (A/R  =  2.0) 


Fig.  III-31  Vortex  shedding  patterns  at  different  aspect  ratios  (Re  =  1. 0  x  104) 


Time  averaged  drag  force  coefficient  is  about  1.07  to  1.15  and  increases  a  little  bit  as  Reynolds 
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number  increases  as  shown  in  Fig.  III-32.  To  be  expected,  time  averaged  lift  or  side  force  coefficient  is 
almost  zero.  Finally,  it  is  observed  that  the  time  averaged  force  coefficients  do  not  strongly  depend  on  the 
aspect  ratio  of  the  geometry  or  the  Reynolds  number. 


Fig.  III-32 
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Time  averaged  force  coefficient  versus  Reynolds  number  at  different  aspect  ratios 


The  plots  in  Figs.  Ill — 33  and  III— 34  demonstrate  the  time  averaged  streamwise  velocity  component 
profiles  at  several  different  distances  from  the  body.  The  time  averaged  streamwise  velocity  component  is 
normalized  based  on  the  free-stream  velocity  and  plotted  along  the  positive  z-direction  at  the  middle 
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section  of  the  height,  i.e.  x  —  z  plane.  As  we  expected,  the  wake  region  behind  a  rectangular 
parallelepiped  is  larger  when  the  aspect  ratio  is  larger  as  shown  in  Fig.  Ill— 33.  Different  Reynolds 
numbers  make  no  difference  in  the  velocity  profiles  as  shown  in  Fig.  III-34.  Namely,  the  wake  region 
behind  a  rectangular  parallelepiped  is  hardly  affected  by  the  Reynolds  number. 
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Fig.  III-34  Velocity  profiles  of  time  averaged  streamwise  velocity  component  at  A/R  =  2.0 


3.  Summary 

The  formation  mechanism  of  the  hairpin-shaped  vortex  shedding  and  the  unsteady  oscillating  flow 
phenomena  are  numerically  investigated. 

In  general,  the  location  of  the  vortex  shedding  is  changed  from  time  to  time  and  rotates  slowly  and 
irregularly  along  the  circumference  of  the  body.  This  irregular  rotating  manner  strongly  depends  on  the 
aspect  ratio  and  Reynolds  number,  and  becomes  weaker  as  the  aspect  ratio  increases  or  as  the  Reynolds 
number  increases  even  though  it  becomes  stronger  as  the  Reynolds  number  increases  when  the  aspect 
ratio  is  1.0. 

The  drag  force  coefficient  has  the  highest  frequency  and  the  least  amplitude.  Frequency  of  the  drag 
force  coefficient  strongly  depend  on  the  aspect  ratio  or  Reynolds  number,  whereas  the  amplitude  of  the 
drag  force  coefficient  does  not  strongly  depend  on  the  aspect  ratio  or  Reynolds  number. 

The  lift  force  coefficient  has  the  largest  amplitude  no  matter  what  the  aspect  ratio  and  the  Reynolds 
number  are.  The  amplitude  of  the  lift  force  coefficient  strongly  depends  on  the  aspect  ratio  and  Reynolds 
number.  The  frequency  of  the  lift  force  coefficient,  however,  does  not  strongly  depend  on  the  aspect  ratio 
or  Reynolds  number. 

The  amplitude  of  the  side  force  coefficient  becomes  weaker  as  the  aspect  ratio  increases  or  as  the 
Reynolds  number  increases.  The  frequency  of  the  side  force  coefficient  is  less  than  that  of  the  lift  force 
coefficient  no  matter  what  the  aspect  ratio  and  the  Reynolds  number  are,  and  does  not  strongly  depend  on 
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the  aspect  ratio  or  Reynolds  number. 

The  size  of  wake  region  depends  on  the  aspect  ratio  of  the  rectangular  parallelepiped.  However,  the 
Reynolds  number  does  not  affect  the  size  of  wake  region. 

Finally,  we  can  conclude  that  the  aspect  ratio  and  Reynolds  number  do  not  strongly  affect  the 
oscillating  lift  or  side  force  coefficient  frequency  although  they  strongly  affect  the  drag  force  coefficient 
frequency.  In  addition,  they  strongly  affect  the  amplitude  of  oscillating  lift  and  side  force  coefficient 
although  they  do  not  strongly  affect  the  amplitude  of  drag  force  coefficient.  They  also  affect  the  vortex 
shedding  location,  which  rotates  slowly  and  irregularly  about  an  axis  through  the  center  of  the  rectangular 
parallelepiped.  However,  they  do  not  affect  the  time  averaged  force  coefficients  such  as  drag,  lift,  and  side 
forces.  The  wake  region  is  hardly  affected  by  the  Reynolds  number. 


III-3.  Analysis  of  Flowfield  Past  a  Circular  Cylinder 
1.  Vortex  shedding  formation 

The  vortex  shedding  patterns  behind  a  circular  cylinder  strongly  depend  on  the  aspect  ratio  of  the 
geometry.  Here,  the  aspect  ratio  of  a  circular  cylinder  is  defined  as  the  ratio  of  width(W)  to  diameter(D) 
of  a  circular  cylinder.  Two  different  vortex  shedding  patterns  are  observed;  one  is  the  symmetric-pair 
vortex  shedding  pattern  shown  in  Fig.  Ill— 35  and  the  other  is  the  alternating  vortex  shedding  pattern 
shown  in  Fig.  Ill — 36. 


(a)Right  side  view  (Re  =  1.0xl03) 


(c)  Right  side  view  ( Re  =  1.0  x  104) 


(b)  Top  view  ( Re  =  1.0  x  103) 


(d)  Top  view  (Re  =  1.0  x  104) 


Fig.  III-35  Vortex  shedding  pattern  at  A/R  =  1. 0 
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(a)  Right  side  view  (Re  —  1.0  x  103) 


(c)  Right  side  view  (Re  —  1.0  x  104) 


(b)  Top  view  (Re  —  1.0  x  103)  (d)  Top  view  (Re  —  1.0  x  104) 

Fig.  III-36  Vortex  shedding  pattern  at  A/R  —  2.0 


(b)  Hairpin-shaped  vortex  (A/R  —  2.0)  (d)  Velocity  distributions  (A/R  —  2.0) 

Fig.  III-37  Flow  features  near  the  body  (right  side  view  at  middle  of  section,  Re  =  1. 0  x  104) 


-49- 


(b)  Hairpin-shaped  vortex  ( A/R  =  2.0)  (d)  Velocity  distributions  (A/R  =  2.0) 

Fig.  III-38  Flow  features  near  the  body  (top  view  at  middle  of  section,  Re  =  1. 0  x  104) 


When  the  aspect  ratio  is  1.0,  the  vortex  shedding  at  the  top  and  bottom  surface  occurs  almost 
symmetrically  with  respect  to  the  horizontal  plane  i.e.  x  —  z  plane,  meanwhile  the  vortex  shedding 
occurs  alternately  at  the  top  and  bottom  surface  when  the  aspect  ratio  is  2.0.  Figures  III— 37  and  III— 3 8 
show  the  vortex  distributions  and  velocity  distributions  after  the  top  vortex  sheds  from  the  body.  The 
upper  part  and  lower  part  of  the  wake  region  are  asymmetric  velocity  distributions  with  respect  to  the 
horizontal  x  —  z  plane  as  shown  in  Fig.  Ill— 37(d)  when  the  aspect  ratio  is  2.0.  Namely,  the  upper  part  of 
the  flow  region  flows  faster  to  the  downstream  direction  than  the  lower  part  of  the  flow  region  so  that  the 
upper  part  of  the  vortex  loop,  which  has  detached  from  the  body  surface,  flows  faster  to  the  downstream 
direction  in  spreading  out.  On  the  other  hand,  the  lower  part  of  the  vortex  loop  flows  to  the  downstream 
direction  slowly  by  forming  the  two  long  legs  shape  of  vortex  because  of  inward  flow  direction  as  shown 
in  Fig.  Ill-  38(d).  When  the  aspect  ratio  is  1.0,  the  upper  part  and  lower  part  of  the  wake  region  have 
almost  symmetric  velocity  distributions  with  respect  to  the  horizontal  x  —  z  plane  as  shown  in  Fig. 
III-37(c)  so  that  the  upper  and  lower  part  of  the  vortex  loop  flow  to  the  downstream  direction  at  the 
almost  same  velocity  magnitude,  which  means  a  symmetric -pair  vortex  shedding  pattern  has  occurred 
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with  respect  to  the  horizontal  plane. 

The  reason  why  there  exist  two  different  vortex  shedding  patterns  is  that  the  side  flow  effect.  Here, 
the  side  flow  effect  means  the  geometrically  three-dimensional  effect  and  becomes  weaker  as  the  aspect 
ratio  increases.  In  other  words,  the  alternating  vortex  shedding  at  the  top  and  bottom  surface  occurs  more 
easily  and  stably  when  the  most  fluid  flow  to  the  downstream  direction  instead  of  flowing  to  the  z- 
direction. 
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(a)  Velocity  distributions  of  u-component  (right  side  view) 
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(b)  Velocity  distributions  of  w-component  (right  side  view) 
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(c)  Velocity  distributions  of  w-component  (left  side  view) 
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(d)  Velocity  distributions  of  w-component  (slightly  tilted 
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right  side  view) 
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(e)  Velocity  distributions  of  w-component  (slightly  tilted  left  side  view) 

Fig.  III-39  Velocity  component  distributions  on  the  vortex  structure  (A/R  =  1. 0,  Re  =  1. 0  x  104) 


Figures  III-39  explains  the  mechanism  of  vortex  shedding  formation  in  a  different  way  using  Q- 
criterion  with  the  velocity  component  distributions  on  the  vortex  structures  when  the  aspect  ratio  is  1.0 
and  Reynolds  number  is  1.0  X  104.  Here,  we  can  understand  the  mechanism  of  symmetric -pair  vortex 
shedding  formation  more  clearly.  The  vortex  loop  simultaneously  detached  from  the  top  and  bottom 
surfaces  sheds  to  the  downstream  direction  as  a  hairpin-shaped  vortex  pair.  At  the  same  time,  the  vortex 
tubes  located  near  the  x  —  z  horizontal  plane  approach  to  each  other  because  of  the  inward  directional 
velocity  distributions  as  shown  in  Fig.  III-39(b)&(c).  Finally,  two  symmetric  vortex  loops  are  formed  and 
flow  to  the  downstream  direction  independently  as  shown  in  Fig.  III-39(a).  Two  vortex  loops  become 
asymmetry  as  they  flow  to  the  downstream. 

Figure  III-40  explains  the  mechanism  of  the  vortex  shedding  formation  at  the  aspect  ratio  of  2.0  at 
Re  =  1.0  x  104.  The  upper  part  of  the  vortex  loop  (colored  by  orange)  flows  faster  to  the  downstream 
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direction  than  the  lower  part  of  the  vortex  loop  (colored  by  green)  as  shown  in  Fig.  III-40(a).  Figure  III- 
40(b)  explains  the  reason  for  the  decreasing  of  the  distance  between  two  long  legs  by  w-velocity 
component  distributions. 
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(a)  Velocity  distributions  of  u-component  (right  side  view) 
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(b)  Velocity  distributions  of  w-component  (bottom  view) 

Fig.  m-40  Velocity  component  distributions  on  the  vortex  structure  (A/R  —  2.0,  Re  —  l.Ox  104) 


(a)Time  history  of  lift  force  coefficient 


(b)  Oil  flow  pattern  at  T*  =  100 
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(c)  Oil  flow  pattern  at  T*  —  160  (d)  Oil  flow  pattern  at  T*  —  200 

Fig.  III-41  Oil  flow  patterns  at  different  time  steps  (A/R=1.0,  Re  =  1. 0  x  103) 


(a)Time  history  of  lift  force  coefficient 


(b)  Oil  flow  pattern  at  T*  —  140 


(d)  Oil  flow  pattern  at  T*  —  180 

steps  ( A/R=2.0 ,  Re  =  1. 0  x  103) 


(c)  Oil  flow  pattern  at  T*  —  160 

Fig.  III-42  Oil  flow  patterns  at 


different 


time 


Figures  III— 4 1  and  III— 42  show  oil  flow  patterns  on  the  body  surface  at  different  time  steps.  When 
the  vortex  shedding  occurs  symmetrically  with  respect  to  the  horizontal  plane,  the  oil  flow  pattern  on  the 
side  wall  surface  also  shows  symmetric  pattern  even  though  the  time  elapse  as  shown  in  Fig.  Ill— 41. 

When  the  vortex  shedding  occurs  alternately  at  the  top  and  bottom  surface,  the  oil  flow  pattern  also 
shows  asymmetric  pattern  as  time  elapses  as  shown  in  Fig.  III-42. 
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In  general,  the  location  of  the  vortex  shedding  rotates  irregularly  about  an  axis  parallel  to  the 
downstream  direction  through  the  center  of  the  body.  However,  no  rotating  tendency  is  observed  for  the 
flow  past  a  circular  cylinder  no  matter  what  the  aspect  ratios  or  Reynolds  numbers  are. 


Fig.  III-43  Time  history  of  the  unsteady  force  coefficient  (A/R  =1.0,  Re  =  1. 0  x  104) 


(c)  T*  =  90 


(d)  T*  =  110 
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Fig.  III-44  Pressure  contours  at  the  location  of  x/D  —  0. 6  (A/R  =1.0,  Re  —  1. 0  x  104) 
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Figure  III— 44  plots  the  pressure  contours  on  the  y  —  z  plane,  where  it  locates  at  x/D  =  0.6  from 
the  center  of  the  body,  at  different  time  shown  in  Fig.  Ill— 43.  The  pressures  are  distributed  evenly  on  the 
plane  at  all  the  different  time  and  are  distributed  symmetrically  with  respect  to  the  horizontal  x  —  z  plane. 
Therefore,  we  can  say  that  the  location  of  the  vortex  shedding  does  not  rotate  about  an  axis  parallel  to  the 
downstream  through  the  center  of  the  body  and  the  vortex  shedding  occurs  at  the  top  and  bottom  surface 
simultaneously  as  a  pair  vortex. 


Fig.  III-45  Time  history  of  unsteady  force  coefficient  (A/R= 2.0,  Re  =  1. 0  x  104) 
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Fig.  III-46 
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Pressure  contours  at  the  location  of  x/D  =  0. 6  (A/R  =2.0,  Re  =  1. 0  x  104) 


On  the  other  hand,  the  pressure  distributions  on  the  y  —  z  plane  at  the  aspect  ratio  of  2.0  are 
symmetry  with  respect  to  the  vertical  x  —  y  plane  at  each  different  time  shown  in  Fig.  Ill— 45  so  that  there 
is  no  symptom  which  demonstrates  the  rotation  of  vortex  shedding  location.  Therefore,  we  can  say  that 
there  is  no  rotating  tendency  of  vortex  shedding  location  and  the  vortex  shedding  occurs  alternately  from 
the  top  and  bottom  surface. 

The  flow  pattern  near  a  circular  cylinder  is  somewhat  different  with  the  flow  pattern  near  a  sphere  or 
near  a  rectangular  parallelepiped.  The  flow  pattern  near  a  circular  cylinder  has  a  flow  pattern  near  a 
curved  surface  plus  a  flow  pattern  near  a  flat  surface.  The  flow  separation  occurs  on  the  curved  surface  as 
well  as  the  front  edge  of  the  side  wall.  Figure  III— 47  shows  the  various  flow  patterns  occurred  near  a 
circular  cylinder. 


(a)  Recirculation  flows  along  the  surface  (b)  Vortex  flows  on  and  near  the  surface  (right  side  view) 


(c)  Vortex  distributions  (front  view)  (d)  Vortex  flows  on  and  near  the  surface  (top  view) 
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(e)  Streamlines  near  the  top  surface 

Fig.  III-47  Flow  features  at  A/R  =  1. 0  and  Re  =  1. 0  x  103 


Figure  III— 48  shows  three  different  types  of  flow  paths  at  a  certain  instantaneous  time;  (1)  Some  of 
inflow  flow  along  the  curved  surface  to  the  downstream  direction  then  recirculate  in  the  wake  region  and 
flow  back  to  the  side  wall,  recirculate  then  flow  to  the  downstream  direction  again  as  shown  in  Fig. 
III-48(a).  (2)  Some  of  inflow  flow  along  the  curved  surface  then  flow  to  the  transverse  direction  and 
recirculate  near  a  side  wall  and  flow  to  the  downstream  direction  as  shown  in  Figs.  Ill— 48(b)  and  (c).  (3) 
Some  of  inflow  bypass  the  curved  surface  and  recirculate  in  the  wake  region  as  shown  in  Fig.  Ill-  48(d). 


(c)  Case  C  (d)  Case  D 

Fig.  III-48  Streamlines  near  the  surface  at  A/R  =  1. 0  and  Re  =  1. 0  x  103 
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Figure  III— 49  plots  the  velocity  distributions  on  the  x-y  plane  at  different  distances  from  the  side 
wall  at  the  Reynolds  number  of  1.0  x  103  and  aspect  ratio  of  1.0.  There  is  no  more  recirculating  zone  at 
z  =  0.60D. 


(c)  Velocity  distributions  at  z  =  0.S7D 


(d)  Velocity  distributions  at  z  =  0.60D 


W 


Fig.  III-49  Velocity  distributions  at  different  locations  (A/R  =  1. 0  and  Re  =  1. 0  x  103) 


2.  Parametric  study 

Figure  III— 50  shows  the  Strouhal  number  variations  according  to  the  variations  of  Reynolds  numbers 
at  different  aspect  ratios.  The  Strouhal  number  of  the  drag  force  coefficient  is  higher  than  that  of  the  lift 
force  coefficient  or  side  force  coefficient  for  all  the  aspect  ratios  and  Reynolds  numbers.  In  additions,  the 
Strouhal  number  of  the  lift  force  coefficient  is  higher  than  that  of  the  side  force  coefficient.  The  Strouhal 
number  of  the  side  force  coefficient  decreases  as  the  aspect  ratio  increases  at  all  the  Reynolds  numbers 
because  the  side  flow  effect  becomes  weaker. 
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The  high  mode  Strouhal  number  for  the  drag  force  coefficient  is  observed  at  the  Reynolds  number  of 
1.0  X  104  for  all  the  aspect  ratios  as  shown  in  Fig.  Ill— 5 1 .  However,  it  is  not  observed  at  the  lower 
Reynolds  number  ranges  including  1.0  X  103, 2.0  x  103,  and  4.0  x  103. 

The  high  mode  Strouhal  number  for  the  lift  and  side  force  coefficient  is  not  observed  at  all  the 
Reynolds  numbers  and  all  the  aspect  ratios. 

In  general,  the  amplitude  of  the  lift  force  coefficient  of  low  frequency  mode  is  much  greater  than  that 
of  the  drag  force  coefficient  of  low  frequency  mode  as  shown  in  Figure  III-52.  In  addition,  the  amplitude 
of  the  drag  force  coefficient  of  low  frequency  mode  is  greater  than  that  of  the  drag  force  coefficient  of 
high  frequency  mode  as  shown  in  Figure  III— 51.  If  there  exists  a  high  frequency  mode  of  the  lift  force 
coefficient,  the  amplitude  of  the  lift  force  coefficient  of  high  frequency  mode,  which  is  directly  related 
with  the  amplitude  of  the  drag  force  coefficient  of  high  frequency  mode  ,  is  much  smaller  than  that  of  the 
lift  force  coefficient  of  low  frequency  mode  so  that  the  high  mode  Strouhal  number  of  the  lift  force 
coefficient  cannot  be  observed  in  FFT. 

Unlike  the  streamwise  directional  flow  or  the  upward  /  downward  directional  flow,  the  transverse 
directional  flow  does  not  include  any  vortex  sheet,  which  is  detached  from  the  curved  body  surface  and 
causes  the  high  mode  Strouhal  number,  so  that  no  high  mode  Strouhal  number  of  the  side  force  coefficient 
is  observed  in  FFT. 

From  now  on,  we  call  the  “low  mode  Strouhal  number”  for  the  drag  force  coefficient  as  simply 
“Strouhal  number”. 
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(c)  A/R  =  2.0 

Fig.  III-51  Fast  Fourier  Transform  of  oscillating  drag  force  coefficient  at  Re  =  1. 0  x  104 

Figure  III— 52  plots  the  amplitude  of  force  coefficient  versus  the  Reynolds  number  at  the  different 
aspect  ratios.  The  amplitude  of  drag  force  coefficient  has  the  minimum  magnitude  among  the  three 
directional  force  coefficients  no  matter  what  the  aspect  ratios  and  Reynolds  numbers  are  as  might  be 
expected.  The  amplitude  of  the  lift  force  coefficient  is  larger  than  that  of  the  side  force  coefficient  at  the 
same  Reynolds  number  when  the  aspect  ratios  are  1.5  and  2.0,  whereas  the  amplitude  of  the  side  force 
coefficient  is  larger  than  that  of  the  lift  force  coefficient  at  the  aspect  ratio  of  1.0  because  the  vortex 
shedding  at  the  top  and  bottom  surface  occurs  almost  symmetrically  with  respect  to  the  horizontal  plane, 
which  means  the  amplitude  of  lift  force  coefficient  becomes  weaker.  As  the  aspect  ratio  increases,  the 
amplitude  of  side  force  coefficient  decreases  at  the  Reynolds  numbers  of  2.0  x  103, 4.0  x  103,  and  1.0  x 
104  because  the  side  flow  effect  becomes  weaker. 
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Fig.  III-52  Amplitude  versus  Reynolds  number  at  different  aspect  ratios 


Fig.  III-53  Time  averaged  force  coefficient  versus  Reynolds  number  at  different  aspect  ratios 

Figure  III— 53  compares  the  time  averaged  force  coefficient.  Time  averaged  drag  force  coefficient  is 
about  0.61  to  0.70  and  increases  as  the  aspect  ratio  increases  at  the  same  Reynolds  number.  On  the  other 
hand,  time  averaged  lift  and  side  force  coefficients  are  almost  zero  values  at  all  the  Reynolds  numbers  and 
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all  the  aspect  ratios. 


Fig.  III-54  Velocity  profiles  of  time  averaged  streamwise  velocity  component  at  Re  =  1. 0  x  103 


Figure  III— 54  plots  the  velocity  profiles  of  time  averaged  streamwise  velocity  component  along  the 
z-direction  at  different  distances  from  the  body  at  the  Reynolds  number  of  1.0  X  103.  The  recirculation 
area  on  the  x  —  z  plane  becomes  larger  as  the  aspect  ratio  increases. 
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Fig.  III-55  Velocity  profiles  of  time  averaged  streamwise  velocity  component  at  A/R  =  2. 0 


Figure  III— 55  plots  the  velocity  profiles  of  time  averaged  streamwise  velocity  component  along  the 
z-direction  at  different  distances  from  the  body  at  the  aspect  ratio  of  2.0.  The  velocity  profiles  hardly 
depend  on  the  Reynolds  numbers. 

Figure  III— 56  plots  the  velocity  profiles  of  time  averaged  streamwise  velocity  component  along  the 
y-direction  at  different  distances  from  the  body  at  the  Reynolds  number  of  1.0  X  103.  The  recirculation 
region  on  the  x  —  y  plane  disappears  faster  as  the  aspect  ratio  decreases. 
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Fig.  III-56  Velocity  profiles  of  time  averaged  streamwise  velocity  component  at  Re  =  1. 0  x  103 


Figure  III— 57  plots  the  velocity  profiles  of  time  averaged  streamwise  velocity  component  along  the 
y-direction  at  different  distances  from  the  body  at  the  aspect  ratio  of  2.0.  The  velocity  profiles  almost  does 
not  depend  on  the  Reynolds  numbers. 
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Fig.  III-57  Velocity  profiles  of  time  averaged  streamwise  velocity  component  at  A/R  =  2. 0 


3.  Summary 

The  formation  mechanism  of  the  hairpin-shaped  vortex  shedding  and  the  unsteady  oscillating  flow 
phenomena  are  numerically  investigated. 

In  general,  the  location  of  the  vortex  shedding  is  changed  from  time  to  time  and  rotates  slowly  and 
irregularly  about  an  axis  through  the  center  of  the  body.  However,  no  rotating  manner  is  observed  for  the 
unsteady  flow  past  a  circular  cylinder. 

Unlike  the  vortex  shedding  behind  a  sphere  or  a  rectangular  parallelepiped,  two  different  vortex 
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shedding  patterns  are  observed  depending  on  the  aspect  ratio  of  geometry;  one  is  the  symmetric-pair 
vortex  shedding  pattern  at  the  aspect  ratio  of  1 .0  and  the  other  is  the  alternating  vortex  shedding  pattern  at 
the  aspect  ratio  of  1.5  and  2.0  no  matter  what  the  Reynolds  number  is. 

The  drag  force  coefficient  has  high  mode  Strouhal  number  as  well  as  low  mode  Strouhal  number. 
The  high  mode  Strouhal  number  is  observed  only  when  the  Reynolds  number  is  1.0  X  104  for  all  the 
aspect  ratios.  The  (low  mode)  Strouhal  number  of  the  drag  force  coefficient  is  higher  than  that  of  lift  force 
coefficient  or  side  force  coefficient  for  all  the  aspect  ratios  and  Reynolds  numbers.  The  amplitude  of  the 
drag  force  coefficient  has  the  minimum  magnitude  among  the  three  directional  force  coefficients  no 
matter  what  the  aspect  ratio  and  Reynolds  number  are.  Time  averaged  drag  force  coefficient  is  about  0.61 
to  0.70  and  increases  as  the  aspect  ratio  increases  at  the  same  Reynolds  number. 

The  lift  force  coefficient  has  only  low  mode  Strouhal  number.  The  Strouhal  number  of  the  lift  force 
coefficient  is  higher  than  that  of  side  force  coefficient.  The  amplitude  of  lift  force  coefficient  is  larger  than 
that  of  the  side  force  coefficient  at  the  same  Reynolds  numbers  when  the  aspect  ratios  are  1.5  and  2.0, 
whereas  the  amplitude  of  lift  force  coefficient  is  smaller  than  that  of  the  side  force  coefficient  at  the 
aspect  ratio  of  1.0.  Time  averaged  lift  force  coefficient  is  almost  zero  value  at  all  the  Reynolds  numbers 
and  all  the  aspect  ratios. 

The  side  force  coefficient  has  only  low  mode  Strouhal  number.  The  Strouhal  number  of  the  side 
force  coefficient  is  the  lowest  no  matter  what  the  aspect  ratios  and  the  Reynolds  numbers  are.  The 
Strouhal  number  of  the  side  force  coefficient  decreases  as  the  aspect  ratio  increases  at  all  the  Reynolds 
numbers.  The  amplitude  of  the  side  force  coefficient  is  smaller  than  that  of  the  lift  force  coefficient  except 
at  the  aspect  ratio  of  1.0  for  all  the  Reynolds  numbers.  As  the  aspect  ratio  increases,  the  amplitude  of  side 
force  coefficient  decreases  at  the  Reynolds  numbers  of  2.0  x  103,4.0  x  103,and  1.0  x  104.  Time 
averaged  side  force  coefficient  is  almost  zero  value  at  all  the  Reynolds  numbers  and  all  the  aspect  ratios. 

The  size  of  recirculation  region  becomes  larger  as  the  aspect  ratio  increases.  In  addition,  it  disappears 
faster  along  the  downstream  direction  as  the  aspect  ratio  decreases.  The  Reynolds  number  does  not  affect 
the  size  of  recirculation  region. 

Finally,  we  can  conclude  that  there  are  two  different  vortex  shedding  patterns  depending  on  the 
aspect  ratio  of  the  geometry  and  there  is  no  rotating  manner  in  the  vortex  shedding.  The  aspect  ratio  of 
geometry  affects  the  amplitude  and  Strouhal  number  of  the  side  force  coefficient,  and  the  size  of  wake 
region.  The  Reynolds  number  affects  the  existence  of  the  high  mode  Strouhal  number  of  the  drag  force 
coefficient. 

III-4.  Overall  Analysis 

In  this  section,  the  Strouhal  number  and  the  amplitude  of  the  oscillating  forces,  and  time  averaged 
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drag  force  coefficient  are  compared  for  the  three  different  geometries  such  as  sphere,  rectangular 
parallelepiped,  and  circular  cylinder. 

The  shape  of  a  sphere  is  approximately  similar  to  the  shape  of  a  circular  cylinder  whose  aspect  ratio 
is  1.0  so  that  the  tendency  of  the  numerical  results  is  similar  to  each  other  as  shown  in  Fig.  Ill— 5  8 . 

A  sphere  has  the  highest  Strouhal  number  for  the  drag  force,  lift  force,  and  side  force  coefficients. 
The  Strouhal  number  of  the  lift  force  coefficient  for  a  circular  cylinder  whose  aspect  ratio  is  1.0  is  as  high 
as  the  Strouhal  number  of  the  lift  force  coefficient  for  a  sphere. 
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Fig.  III-58  Strouhal  number  versus  geometry 


A  rectangular  parallelepiped  has  the  lowest  Strouhal  number  for  the  lift  force  and  side  force 
coefficients.  Namely,  a  sphere  is  the  geometric  shape  which  oscillates  most  frequently  and  a  rectangular 
parallelepiped  is  the  geometric  shape  which  oscillates  most  infrequently. 
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Fig.  III-59  Amplitude  versus  geometry 


Figure  III— 59  compares  the  amplitude  versus  geometry.  The  amplitude  of  the  drag  force  coefficient 
is  generally  very  small  so  that  it  may  be  useless  to  investigate  which  geometry  has  the  largest  amplitude 
even  though  a  rectangular  parallelepiped  has  the  largest  amplitude. 

A  rectangular  parallelepiped  has  the  largest  amplitude  of  the  lift  force  coefficient.  The  amplitude  of 
the  side  force  coefficient  decreases  as  the  aspect  ratio  increases.  If  we  confine  the  case  which  has  the 
aspect  ratio  of  1 .0,  a  rectangular  parallelepiped  has  the  largest  amplitude  of  the  side  force  coefficient.  And, 
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a  sphere  has  the  smallest  amplitude  of  drag  force,  lift  force,  and  side  force  coefficient. 
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Fig.  III-60  Time  averaged  drag  force  coefficient  versus  different  geometry 

It  is  very  clear  that  a  rectangular  parallelepiped  has  the  largest  time  averaged  drag  force  coefficient 
and  a  sphere  has  the  smallest  drag  force  coefficient  as  might  be  expected  as  shown  in  Fig.  Ill— 60. 
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IV.  CONCLUSIONS 

The  unsteady  flowfield  past  a  simplified  MAV  shaped  geometry  such  as  sphere,  rectangular 
parallelepiped,  and  circular  cylinder  has  been  analyzed  numerically  using  Fluent  to  prepare  a  design 
guideline  for  MAV.  The  mechanism  of  vortex  shedding,  the  Strouhal  number  and  the  amplitude  of  the 
oscillating  drag,  lift,  and  side  force  coefficient  are  studied  in  the  range  of  Reynolds  numbers  from 
1.0  X  103  to  1.0  X  104  with  three  different  aspect  ratios  of  1.0,  1.5,  and  2.0. 

Even  though  the  mechanism  of  vortex  shedding  formation  is  very  similar  to  each  other,  the  vortex 
shedding  pattern  has  some  differences.  There  exists  a  high  mode  Strouhal  number  of  the  drag  force 
coefficient  for  the  flow  past  a  sphere.  However,  there  is  no  high  mode  Strouhal  number  for  the  flow  past  a 
rectangular  parallelepiped.  In  case  of  flow  past  a  circular  cylinder,  there  is  a  high  mode  Strouhal  number 
for  the  drag  force  coefficient  only  at  the  Reynolds  number  of  1.0  X  104.  Actually,  the  amplitude  of  high 
frequency  mode  is  very  tiny  compared  with  other  force  coefficient  amplitudes  even  though  it  has  very 
highly  oscillating  frequency. 

In  general,  the  vortex  sheds  form  the  body  in  randomly  rotating  pattern  about  an  axis  parallel  to  the 
downstream  direction  through  the  center  of  the  body.  However,  no  rotating  manner  is  observed  for  the 
flow  past  a  circular  cylinder. 

The  vortex  sheds  alternately  from  the  top  and  bottom  surface  for  the  flow  past  a  sphere  and  a 
rectangular  parallelepiped.  It  also  sheds  alternately  from  the  top  and  bottom  surface  for  the  flow  past  a 
circular  cylinder  when  the  aspect  ratios  are  1.5  and  2.0.  However,  a  symmetric-pair  vortex  shedding 
occurs  at  the  flow  past  a  circular  cylinder  when  the  aspect  ratio  is  1 .0  in  the  range  of  Reynolds  numbers 
from  1.0  x  103  to  1.0  x  104. 

Generally,  the  drag  force  coefficient  has  the  highest  Strouhal  number  meanwhile  it  has  the  least 
amplitude.  The  amplitude  of  the  lift  force  coefficient  is  the  largest  except  for  the  flow  past  a  circular 
cylinder  at  the  aspect  ratio  of  1.0.  The  amplitude  of  the  side  force  coefficient  decrease  as  the  aspect  ratio 
increases  for  the  flows  past  a  rectangular  parallelepiped  and  a  circular  cylinder.  The  Strouhal  number  of 
the  side  force  coefficient  is  the  lowest  in  general. 

A  sphere  has  the  highest  Strouhal  number  for  the  drag  force,  lift  force,  and  side  force  coefficients 
where  as  a  rectangular  parallelepiped  has  the  lowest  Strouhal  number  for  the  drag  force,  lift  force,  and 
side  force  coefficients.  A  rectangular  parallelepiped  has  the  largest  amplitudes  of  the  drag  force,  lift  force, 
and  side  force  coefficient  whereas  a  sphere  has  the  smallest  amplitudes  of  the  drag  force,  lift  force,  and 
side  force  coefficient.  A  rectangular  parallelepiped  also  has  the  largest  drag  force  coefficient. 

As  an  additional  remark,  a  very  careful  investigation  is  recommended  for  the  flowfield  analysis  past 
a  geometry  such  as  a  rectangular  parallelepiped  or  a  circular  cylinder  whose  aspect  ratio  is  1.0  because  it 
has  strong  three-dimensional  geometric  characters. 
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APPENDIX  A -Sphere 


Fig.  A1  Time  history  of  unsteady  drag  force  coefficient 


(c)  Re  =  4,000  (d)  Re  =  10,000 

Fig.  A2  Power  spectrum  of  drag  force  coefficient 
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Fig.  A3  Time  history  of  unsteady  lift  force  coefficient 
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Fig.  A4  Power  spectrum  of  lift  force  coefficient 
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Fig.  A5  Time  history  of  unsteady  side  force  coefficient 
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Fig.  A6  Power  spectrum  of  side  force  coefficient 
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Fig.  A7  Oil  flows  on  the  top  and  bottom  surface  with  same  legend  of  skin  friction  coefficient 
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Fig.A8  Vortex  structure  (side  view  and  top  view) 
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APPENDIX  B  -  Rectangular  Parallelepiped 


(a)  Drag  force  coefficient 


(b)  Lift  force  coefficient 


(c)  Side  force  coefficient 


Fig.  B1  Time  history  of  unsteady  force  coefficient  at  Re  =  1, 000 
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(a)  Drag  force  coefficient  (b)  A/R  —  1.0 


(c)  A/R  =  1.5  (d)  A/R  =  2.0 


Fig.  B2  Power  spectrum  of  drag  force  coefficient  at  Re  —  1, 000 
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(a)  Lift  force  coefficient  (b)  A/R  =  1.0 
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Fig.  B3  Power  spectrum  of  lift  force  coefficient  at  Re  —  1, 000 
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(a)  Side  force  coefficient  (b)  A/R  —  1.0 


(c)  A/R  =  1.5  (d)  A/R  =  2.0 


Fig.  B4  Power  spectrum  of  side  force  coefficient  at  Re  —  1, 000 
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(a)  A/R  =  1.0 


(c)  A/R  =  2.0 


Fig.  B5  Vortex  structure  (side  view  and  top  view)  at  Re  —  1, 000 
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Fig.  B6  Time  history  of  unsteady  force  coefficient  at  Re  =  2, 000 
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Fig.  B7  Power  spectrum  of  drag  force  coefficient  at  Re  =  2, 000 
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(a)  Lift  force  coefficient  (b)  A/R  —  1.0 
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Fig.  B8  Power  spectrum  of  lift  force  coefficient  at  Re  —  2, 000 
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(a)  Side  force  coefficient  (b)  A/R  —  1.0 
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Fig.  B9  Power  spectrum  of  side  force  coefficient  at  Re  =  2, 000 
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(a)  A/R  =  1.0 


(c)  A/R  =  2.0 


Fig.  BIO  Vortex  structure  (side  view  and  top  view)  at  Re  —  2, 000 
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(a)  Drag  force  coefficient  (b)  A/R  —  1.0 


(c)  A/R  =  1.5  (d)  A/R  =  2.0 


Fig.  B12  Power  spectrum  of  drag  force  coefficient  at  Re  —  4, 000 
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(a)  Lift  force  coefficient  (b)  A/R  —  1.0 


(c)  A/R  =  1.5  (d)  A/R  =  2.0 


Fig.  B13  Power  spectrum  of  Lift  force  coefficient  at  Re  —  4, 000 


-92- 


(a)  Side  force  coefficient  (b)  A/R  —  1.0 


(c)  A/R  =  1.5  (d)  A/R  =  2.0 


Fig.  B14  Power  spectrum  of  Side  force  coefficient  at  Re  —  4, 000 
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Fig.  B15  Vortex  structure  (side  view  and  top  view)  at  Re  =  4, 000 


-94- 


1.30 


1.25 
1.20 
1.15 

1.10 

1.05 

1.00 

0  10  20  30  40 

Time(T*) 

(a)  Drag  force  coefficient 


A/R  =  1 .0 
A/R  =  1 .5 
A/R  =  2.0 


(b)  Lift  force  coefficient 


(c)  Side  force  coefficient 

Fig.  B16  Time  history  of  unsteady  force  coefficient  at  Re  =  10, 000 


-95- 


(a)  Drag  force  coefficient  (b)  A/R  —  1.0 
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Fig.  B17  Power  spectrum  of  drag  force  coefficient  at  Re  —  10, 000 
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(a)  Side  force  coefficient  (b)  A/R  —  1.0 
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Fig.  B18  Power  spectrum  of  lift  force  coefficient  at  Re  —  10, 000 
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(a)  Side  force  coefficient  (b)  A/R  —  1.0 
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Power  spectrum  of  side  force  coefficient  at  Re  —  10, 000 
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Fig.  B20  Vortex  structure  (side  view  and  top  view)  at  Re  =  10, 000 
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APPENDIX  C  -  Circular  Cylinder 
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Fig.  C2  Power  spectrum  of  drag  force  coefficient  at  Re  =  1, 000 
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Fig.  C3  Power  spectrum  of  lift  force  coefficient  at  Re  —  1, 000 
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Fig.  C4  Power  spectrum  of  side  force  coefficient  at  Re  =  1, 000 
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Fig.  C5  Vortex  structure  (side  view  and  top  view)  at  Re  =  1, 000 
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-  106- 


Amplitude 


(a)  Lift  force  coefficient  (b)  A/R  =  1.0 


(c)  A/R  =  1.5  (d)  A/R  =  2.0 

Fig.  C8  Power  spectrum  of  lift  force  coefficient  at  Re  =  2, 000 


-  107  - 


Amplitude 


(a)  Side  force  coefficient  (b)  A/R  =  1.0 


(c)  A/R  =  1.5  (d)  A/R  —  2.0 
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Fig.  CIO  Vortex  structure  (side  view  and  top  view)  at  Re  =  2, 000 
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